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ABSTRACT 

We present single-Sersic two-dimensional model fits to 167, 600 galaxies modelled inde- 
pendently in the ugrizYJHK bandpasses using reprocessed Sloan Digital Sky Survey 
Data Release Seven (SDSS DR7) and UKIRT Infrared Deep Sky Survey Large Area 
Survey (UKIDSS-LAS) imaging data available from the GAMA database. In order to 
facilitate this study we developed SIGMA, an R wrapper around several contemporary 
astronomy software packages including Source Extractor, PSF Extractor and GALFIT 
3. SIGMA produces realistic 2D model fits to galaxies, employing automatic adaptive 
background subtraction and empirical PSF measurements on the fiy for each galaxy 
in GAMA. Using these results, we define a common coverage area across the three 
GAMA regions containing 138, 269 galaxies. We provide Sersic magnitudes truncated 
at 10 Te which show good agreement with SDSS Petrosian and GAMA photometry 
for low Sersic index systems (n < 4), and much improved photometry for high Sersic 
index systems (n > 4), recovering as much as Am = 0.5 magnitudes in the r band. We 
employ a K band Sersic index/w — r colour relation to delineate the massive (n >~ 2) 
early-type galaxies (ETGs) from the late-type galaxies (LTGs) . The mean Sersic index 
of these ETGs shows a smooth variation with wavelength, increasing by 30% from g 
through K. LTGs exhibit a more extreme change in Sersic index, increasing by 52% 
across the same range. In addition, ETGs and LTGs exhibit a 38% and 25% decrease 
respectively in half-light radius from g through K. These trends are shown to arise due 
to the effects of dust attenuation and stellar population/metallicity gradients within 
galaxy populations. 

Key words: galaxies: structure - galaxies: fundamental parameters - astronomical 
data bases: miscellaneous - catalogues 
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1 INTRODUCTION 

The shapes and sizes of galaxies are not random but are 
defined by the orbital motions of their constituent stellar 
populations, arranging themselves into elliptical, bulge, disk, 
and bar like structures. Exactly why and how these struc- 
tures come about is somewhat a mystery which no doubt 
relates to a complex formation history involving collapse, 
merging, infall, secular evolution, and feedback processes as 
well as the precise nature of the coupling between the dark 
matter, gas, dust and stars and the influence of the lar- 
ger halo in which the galaxy might reside (group, cluster, 
etc.), and the broader environment (filament, void, nexus, 
etc.). The combination of variations in, for example, galaxy 
structure, formation history, evolution and relative environ- 
ment lead to distinct measurable effects on global galaxy 
parameters such as colour, concentration and size. The ulti- 
mate goal of structural analysis is to inform this discussion 
by robustly isolating and quantifying these parameters and 
exploring correlations between these properties and those 
obtained by other means, such as dynamical information. 

Once the underlying structure of a galaxy is under- 
stood the overarching morphological class may be determ- 
ined, and from this we can explore correlations with, for 
example, environment through the well known morphology- 
density relation (Dressier 19801, i.e., the apparent prefer- 



ence for red, passive galaxies in the dense cores of galaxy 
groups and clusters. Several mechanisms have been sugges- 
ted to explain this feature, most notably the combined ef- 



1993 



Gott 



fects of strangulation ( Larson et al. 111980 Kauffmann et al 



Diaferio et al.|200T l , ram pressure stripping ( Gunn & 



19721, harassment (Moore et al. 1996 1 and tidal in- 



teractions and merging (Park et al. 20081. Recent studies 



by, e.g., van der Wei (20081; Welikala et al. (2008 20091; 



Bamford et al. ( 2009 1 and others confirm this morphology- 



environment connection; however, they suggest that the re- 
lation between structure and morphology is less apparent. 
Indeed, it appears that the mass of a galaxy rather than 
the environment in which it resides is more influential in de- 
termining its structure, highlighting the importance of stel- 
lar mass estimates. 

As an example of the connection between galaxy struc- 
ture and the physical processes of galaxy formation |Dal- 



canton, Spergel, & Summers (1997 1 and independently Mo, 
Mao, fc WhiteJ ( |1998| , both following on from [Fall fc Ef- 



stathiou (19801, relate the scale length of the disk to the 



angular momentum of a galaxy's dark matter halo. In addi- 
tion, numerous properties of the bulge component are now 
known to relate to the mass of the super-massive black hole 



(e.g.,|Haring fc Rix|2004l [Novak, Faber, fc Dekel|2006[ |Gra 
|ham fc Driver|2007[ ). Variations in structural properties as a 
function of wavelength (e.g.. La Barbera et al.|2010 l enable 
the extraction of colour gradients, potentially implying the 
direction of disk growth (e.g., inside out,|Barden et al.|2005[ 



Bakos et al.|2008[|Trujillo et al. |20"09l | Wang et al.|2011| , or 

arguing for the redistribution of populations from the inner 
to outer regions (;Roskar et al. 2008), possibly coupled with 
bar formation ( jPebattista et al.^2006j . 

The physics underpinning galaxy structure is relatively 
immature, despite the very long history dating back to 
Knox-Shaw, Reynolds and Hubble and essentially consists 
of spot-check simulations which focus on a particular phe- 



nomena in a non-cosmological context (for recent develop- 
ments see iRoskar et al. 2010 and Agertz et al. 2011) . For 



example, numerical models can readily produce bar, pseudo- 
bulge, spiral patterns, and spheroidal structures through 
coupled rotation, secular evolution, shock-wave propagation, 
and merging history. Until recently the very thin nature 
of the spiral disks has presented a particular challenge for 
numerical models, with numerical simulations in particular 
forming small thick disks, mainly because of the high level 
of merging (Navarro fc Steinmetz 19971. This is in stark 



contrast to a number of independent empirical studies (e.g.. 
Driver et al.|2007|[Gadotti 2009l[Tasca fc White|2011[ ) which 



estimate that approximately 60% of the stellar mass in the 
universe today lies within disk systems, suggesting a more 
quiescent merger history (however, see Hopkins 



Quataert 

20 10|on the stability of gas ri ch disks). In addition, studies 



by |Menendez-Delmestre et al. ( |2007[ | suggest that up to 67% 
of spiral galaxies contain a barred structure, further complic- 
ating simulation efforts. However, numerical simulations are 



now starting to produce realistic disk systems (Governato 
etar|2 007 2009 A gertz et aL| 2009 201 1| albeit with heavily 



controlled initial conditions, more quiescent merger histories 
and a greater degree of gas infall. 

Beyond a distance of ~ 100 Mpc detailed structural 
studies have been relatively rare and mostly confined to 
the deep yet very narrow pencil beam surveys from the 
Hubble Space Telescope (HST). It was only following the 
refurbishment of HST that structural analysis once again 
became a mainstream study (Driver et al. 1995a|b 19981. 
HST provides kpc resolution across the full path length of 
the Universe, which is now also becoming possible with AO 
ground-based systems ( Huertas-Company et al!]|2007 |. The 
conjunction of development in numerical models and this 
new ability to resolve the shapes and sizes of galaxies at 
any distance has led to a dramatic renewed interest in struc- 
tural analysis. One interesting claim is the apparent remark- 
able growth of galaxy sizes since intermediate redshifts (e.g., 
Trujillo fc Pohlen|2005l [Harden et al.|2005| [Mcintosh et al] 



2005|[Trujillo et ai.|2006||2007[[Weinzirl et al.|2011| , pote~ 
tially supporting the notion of recent growth in disk systems 
following an earlier aggressive merger phase at z ~ 2 (see 
[Driver et al. [[19961 12005| |. An alternative suggestion which 
does not require galaxy growth through mergers is the trans- 
formation of some of these so-called "red nuggets" into the 
bulges of disk galaxies via the accretion of a cold gas disk 



(Graham 20111 



However, structural analysis is not trivial to implement 
and interpret correctly, and is plagued by a number of key 
issues. In particular: 

(i) Wavelength bias. At different wavelengths, light traces 
varying stellar populations (Block et al. 1999 1. Typically this 



is a young stellar population at shorter wavelengths and 
an older stellar population at longer wavelengths. For this 
reason, it is vital when comparing structural properties to 
compare properties measured at the same rest wavelength. 

(ii) Dust attenuation. Dust is predicted to modify not 
only the recovered total fiux as a function of wavelength 
(e.g., [Tuffs et al^[2004[ [Pierini et aL][2004| but also galaxy 



sizes, shapes and concentrations (sec for example MollenhofF 
et al |2006[ [Graham fc Worley^2008j . Dust can vary enorm- 



ously from system to system with significant environmental 
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dependencies and strong evolution with redshift. Each in- 
dividual galaxy ultimately requires either a dust correction 
or analysis at rest-NIR wavelengths where dust will have a 



smaller impact (see photon escape fraction curve in Driver 



et al. 2008 1 . It is obvious that any attempt to model the dust 



in galaxies raises the larger problem of degeneracies appear- 
ing between additive and subtractive flux components. 

(iii) Local minima during the minimisation process. For a 
single profile fit there are often 7 free parameters, with that 
number rising for multi-component fits. The surface within 
this parameter space is known to be complex, containing 
multiple local minima representing potentially non-physical 
results, e.g., a bulge which contributes significantly more 



flux to the outer regions of a galaxy than the disk (Gra- 
|ham|200l| . Other than manual checks of the output, various 
methods may be employed to reduce the risk of divergence 
on an incorrect result including; constraints applied during 
the minimisation routine and employing an automated lo- 
gical filter (e.g., [Allen et al.|2006[ ). 

(iv) What lies below the limiting isophote. Whilst the 
surface brightness profile of some galaxies behaves as ex- 



pected out to very faint magnitudes (e.g., NGC 300 Bland- 
Hawthorn et"aL]|2005| [Vlajic et al.|[2509l NGC 7793' |Vlajic 
et al. 20ff ), the potential myriad of phenomena present in 
the outer wings of many systems may cause deviations away 
from a typical light profile. These include truncated and 



anti-truncated disks (Erwin et al. 2005 Pohlen & Trujillo 
2006[), UV excesses ( |Bush et al.||2010|l, tidal d ebris, halos 
( [Barker et al ."2009 McConnachie et al.[[2009| an d minor 
merger fossil records (Martinez-Delgado et al.|2010 |. In fact, 
the outer regions of galaxies may defy any systematic pro- 
file fitting into a restricted number of structures. The ac- 
curacy of any estimation of the background sky and gradi- 
ents therein will also no doubt affect analyses of these outer 
structures. 

(v) The number of components required. When consid- 
ering the structure of very nearby galaxies, the deeper one 
looks the more one finds. Some galaxies, even in the dust 



free 3.8 micron bands, require up to six components (Buta 



et al. 2010 1 before a satisfactory fit can be obtained. In many 



cases there is uncertainty as to how many components are 
required, how to quantitatively decide this in an automated 
and repeatable fashion, and which components are funda- 
mental and which perhaps secondary. For example, should 
bar and pseudo-bulge flux be incorporated into a single disk 
model or kept distinct. 

(vi) Sky estimation. Understanding the background sky 
level at the position of your primary object of interest is cru- 
cial in producing meaningful measurements of that galaxy. 
Considerations must be made in regards accuracy and speed 
of estimating the background. 

(vii) Systematic selection bias. Sample bias will be in- 
troduced due to size, resolution, orientation, profile shape 
and smoothing scale limitations ( [Phillipps &: Disney[[1986{ 
Driver|1999 1. Samples of galaxies are usually selected based 
on global criteria, such as magnitude. However, it becomes 
non-trivial to transcribe these global limits into appropri- 
ate limits for galaxy sub-components, e.g., a certain type of 
disk may only have been detected because it also contains a 
prominent bulge. 

There are several publicly-available galaxy modelling codes 



in common usage including GIM2D (Simard et al. 20021, 
BUD DA ( [de Souza et al.|2004|, GASPHO T ( [Pignatelli et al. 
2006[ ) and? GALFIT 3 ( [Peng et al.[[2010| . In addition, there 
are a number of software pipelines, wrappers around contem- 
porary astronomy software, that aim to automate the pro- 
cess of galaxy modelling including GALAPAGOS (Barden 
et al, 2011, submitted) and PyMorph ( [Vikram et al._2010| . 
These packages all have their advantages and disadvantages 
and have been compared in a number of external studies 
(e.g. Haussler et al.[[2007 Hoyos et al.[|20lT l in addition to 
their own internal comparisons, and so we refer the reader 
to these papers for discussions of the pros and cons between 
ID V 2D fitting and the actual minimisation algorithms em- 
ployed. For this body of work GALFIT was chosen for its 
ease of use and high-quality realistic model outputs, plus the 
ability to perform simultaneous modelling of nearby neigh- 
bours to the primary galaxy. 



In this series of papers we introduce and utilise SIGMA, 
an automated code designed to produce single-Sersic and 
multi-component profile fits for galaxies in the GAMA 
dataset. Using SIGMA, this paper presents one of the 
largest catalogues of multi-wavelength single-Sersic model 
fits; 167, 600 galaxies modelled independently across 9 band- 
passes. This catalogue is currently in use to aid in meas- 
urement of the evolution in the size- (stellar mass) distribu- 
tion of galaxies (Baldry et al., 2011); explore star formation 
trends as a function of morphology (Bauer et al., in prep.); 
to further understand the cosmic SED from 0.1 micron to 1 
mm (Driver et al., 2011); to apply dust corrections to galax- 
ies observed at multiple inclinations (Grootes et al., 2012); 
to explore the dust properties and star-formation histories of 



local submillimetre selected galaxies Rowlands et al. (2011 1 



better constrain stellar mass measurements by providing 
total flux corrections ( Taylor et al.[[20TT l; comment on the 
quenching of star formation in the local universe (Taylor et 
al., in prep.); explore the relation between galaxy environ- 
ments and their star formation rate variations (Wijesinghe 
et al., in prep.); provide a new method for automatic mor- 
phological classification (Kelvin et al., in prep.); and further 
explore the relation between environment (i.e., halo mass; 
Robotham et al. 2011 1, morphology and structure (Kelvin 
et al., in prep.). Future studies building on SIGMA will in- 
corporate advanced logical filtering and profile management 
to produce multi-component fits for a low redshift sample, 
allowing full structural decomposition into bulge-disk-bar, 
etc. (Kelvin et al., in prep). 



This paper is organised as follows. We first outline the 
GAMA data in Section [5] We describe the SIGMA (Struc- 
tural Investigation of Galaxies via Model Analysis) pipeline 
developed to process this data and produce robust 2D galaxy 
models in Section [3|and present SIGMA single-Sersic results 
for 167, 600 objects modelled independently in ugrizY J H K 
from the GAMA Phase I study in Section [4] From this 
large catalogue we establish a common coverage sample of 
138, 269 galaxies. Finally, we further explore the wavelength 
dependence on recovered structural parameters in Section[5] 
A standard cosmology of Hq = 70 km s~^ Mpc~^, Q,m = 0.3, 
= 0.7 is assumed throughout. 
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2 DATA 

The GAMA survey is a combined spectroscopic and multi- 
wavelength imaging programme designed to study spatial 
structure in the nearby [z < 0.25) Universe on scales of 
Ikpc to IMpc (see Driver et al.||2009 for an overview). The 
survey, after completion of Phase I, consists of three regions 
of sky each of 4 deg (Dec) x 12 deg (RA), close to the equat- 
orial region, at approximately 9'' (G09), 12'' (G12) and 14.5'' 
(G15) Right Ascension (see Table[l]and FigurejlJ. The three 
regions were selected to enable accurate characterisation of 
the large scale structure over a range of redshifts and with 
regard to practical observing considerations and constraints. 
They lay within areas of sky scheduled for survey by both 
SDSS ( [Abazajian et al.|2009| as part of its Main Survey, and 
UKI RT as part of the UKID SS Large Area Survey (UKIDSS 
LAS; Lawrence et al. 2007 1. These data provide moderate 



depth and resolution imaging data in ugrizY J H K suitable 
for analysis of nearby galaxies. The accompanying spectro- 
scopic input catalogue was derived from the SDSS PHOTO 
param eter [Stoughton et al.| ( |2002| as described in [Baldry 



et al. (20101. The GAMA spectroscopic programme (Ro- 
botham et al.|20T0 1 commenced in 2008 using AAOmega on 
the Anglo- Australian Telescope to obtain distance informa- 
tion (redshifts) for all galaxies brighter than r < 19.8 mag. 
The survey is ~ 99 percent complete to r < 19.4 mag in G09 
and G15 and r < 19.8 mag in G12 (see Table[T] column 6), 
with a median redshift of z ~ 0.2. Full details of the GAMA 
Phase I spectroscopic program, key survey diagnostics, and 
the GAMA public and team databases are given in IDriver] 
et al.|(|2011||. 



The data used in this paper are obtained from the 
GAMA Database (Driver et al. 2011 1, and include repro- 



cessed imaging from the SDSS [ugriz) and UKIDSS LAS 
(YJHK) archive as described in |Hill et alT| ( |20TT| . The re- 
processing involves the creation of large single image mosa- 
ics for each region in each filter, commonly referred to as 
SWARPed images (swpim) due to the SWarp software used 
in their creation ( Berlin et al.|2002 I . Associated weight map 
mosaics {swpwt) are also constructed. The mosaicing process 
is described in full in [Hill et ah] ( |2011| . In brief, all native 
reduced frames are downloaded from the respective archives 
(SDSS DR7 and ROE/WFAU) and scaled to a single uni- 
form zero point. For SDSS, the input data are the correc- 
ted ifpC) Data Release 7 (DR7) files, with the data having 
already been bias subtracted and flat-fielded as part of the 
SDSS frames pipeline (IStought on et al.||2002"| Section 4.4). 
The UKIDSS LAS data has been collected from the UKIDSS 



et al. 


2007a 


et al. 


(2007 



Early Data Release (EDR; Dye et al. 20061 and data re- 
leases 1 and 2 (DRl; [Warren et al.||2007b[ DR2; |Warren 



The UKIDSS project is defined in Lawrence 



UKIDSS uses the United Kingdom Infrared 



Telescope Wide Field Camera (WFCAM; Casali ct al. 2007 1 



The photometric system is described in Hewett et al. ( 2006 1 , 



and the calibration is described in Hodgkin et al. (20091. 
The pipeline processing and science archive are described in 
Irwin et al (in prep.) and Hambly et al. (2008 1. 



Once these input data have been obtained and calib- 
rated, SWarp is then used to combine them into a single 
image mosaic at a resolution of 0.339" arcscconds per pixel 
in the TAN projection system ( jCalabrctta fc Greisen||2002 I 
centred within each GAMA region as appropriate. Note that 



we are using version 2 SWarp mosaics scaled to a slightly 
higher resolution (0.339 " j^th an the version 1 mosaics (0.4") 
described in Hill et al. (2011 1. Version 2 mosaics are a min- 



imum of 193900 X 79700 pixels each, with each individual 
FITS file ^ 60GB in size. The process used to create the 
version 2 mosaics is identical except the regions have been 
expanded in preparation for GAMA Phase II operations and 
at higher resolution in preparation for matching to VISTA 
data in due cours^E] 

As part of the SWarp mosaicing process the back- 
ground is removed on each individual frame prior to merging 
using a 256 x 256 pixel median filtered mesh which itself is 
median filtered within a 3 x 3 mesh. The original SDSS and 
UKIDSS data are typically held in chunks of 2048 x 1489 
pixels and 2072 x 2072 pixels respectively, at comparable 
pixel scales (SDSS: 0.396"/pixel, UKIDSS: 0.4"/pixel). At 
the native pixel resolution the mesh size therefore equates 
to 101.4" X 101.4" and 102.4" x 102.4" respectively and so 
structures with half-light radii less than ~ 17" should be 
unaffected by the background smoothing 

In addition to the science image frames are the asso- 
ciated weight maps. Because of the zero-point normalisa- 
tion across all data, and overlapping edge duplication in 
the SDSS data, the actual weight map values produced by 
SWarp are an approximation of their correct value. How- 
ever, the weight maps remain useful as a record of which 
stars can be associated with which pre-mosaiced frame for 
the purposes of detailed PSF modelling (described in section 
[33. 

To create our sample of galaxies for modelling, we ex- 
tracted 167, 600 galaxies from the GAMA Tiling Catalogue 
version 11 (TilingCatvll), selecting all galaxy-like objects 
using the GAMA catalogue fiag SURVEY_CLASS> iQ 
The output from these galaxies is stored in the catalogue 
SersicCatv7, presented in Section|4] These galaxies, the mo- 
saiced images and the weight maps constitute our input 
dataset and are all available from the GAMA databas^E] 
as: TilingCatvll, x.mosaic.v2.fits and x. weight. v2. fits, where 
ii=ugrizY J H K . 



^ This increased resolution has been chosen to match that which 
is expected for future VISTA VIKING data releases, allowing easy 
cross-wavelength cross-facility comparison of data. Original SDSS 
and UKIDSS resolutions of 0.396" and 0.4" respectively place a 
limit on how high one is able to artificially increase the resolution 
of mosaiced data, requiring increasing amounts of interpolation 
with increasing artificial resolution. Further details may be found 
in Liske et al. (in prep.). 

^ These larger, higher-resolution version 2 mosaics will be 
released shortly via the GAMA website: http;//www.gama- 
survey.com. 

3 UKIDSS J band data and selected UKIDSS EDR fields in both 
H and K bands were microstepped. These data are typically 
stored in chunks of 4103 X 4103 pixels at a native resolution of 
0.2" /pixel, giving a mesh size of 51.2" X 51.2". 

SURVEY_CLASS is a flag present in many GAMA catalogues 
allowing one to quickly select pre-deflned subsets of the GAMA 
data. 

^ The GAMA database can be found at http://www.gama- 
survey.org/ database. 
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Figure 1. Coverage in each bandpass across the three GAMA regions is shown, where white space indicates a lack of coverage at that 
position. Different bandpasses are offset in declination for plotting purposes, with the bottom row (black data points) positioned at the 
correct GAMA coordinates. The five SDSS bands (ugriz) were taken simultaneously and therefore are represented together, whilst the 
four UKIDSS bands (YJHK) are shown separately. The data points are coloured according to the local value of the PSF FWHM at 
that position, giving an indication of the variation in seeing (SDSS data is taken from the r band). See Section |3.4| for further details on 
the derivation of the PSF's. The bottom row represents a common coverage area, containing galaxies that have been observed across all 
nine bands (92.8% of the total). A common coverage sample is derived from this area, and used for all subsequent figures. 



Region 


RA (J2000) 


Dec (J2000) 




Nobj 


^comp 


G09 
G12 

G15 


129°. <a < 141°. 
174°. <a< 186°. 
211°. 5 <a< 223°. 5 


-1°.0 < 5 < 3°.0 
-2°.0 <S < 2°.0 
-2°.0 < (5 < 2°.0 


19.4 (19.8) 
19.8 (19.4) 
19.4 (19.8) 


30, 289 (48, 548) 
50,868 (32,747) 
33,205 (51,217) 


99.23% (79.36%) 
99.12% (99.39%) 
98.95% (79.27%) 



Table 1. GAMA region definitions. The GAMA main survey definitions are based on SDSS extinction corrected r-band Petrosian 
magnitude limits, the depth of which varies between r = 19.4 mag in G09/G15 and r = 19.8 mag in G12. Comparison magnitude limits 
are shown in brackets for reference. Number counts and redshift completeness arc based on objects which passed star-galaxy separation 
in the GAMA TilingCatvll (see |Baldry et al.||20Tol l for further details). 
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3 SIGMA: AUTOMATED GALAXY 
MODELLING 

SIGMA (Structural Investigation of Galaxies via Model 
Analysis) is an automated front-end wrapper which util- 
ises a wide-range of image analysis software and a series 
of logical filters and handlers to perform bulk structural 
analysis on an input catalogue of galaxies. This is primar- 



ily achieved through the use of Source Extractor ( Bertin & 



Arnouts 19961, PSF Extractor (Bertin and Delorme, priv. 
comm.) and GALFIT 3 ( |Peng et al^|2010| |, with additional 
packages also created and utilised to aid in the fitting pro- 
cess. Key to this process is the galaxy modelling software 
GALFIT. GALFIT is able to create a reaUstic model of 
each input galaxy by fitting one or more analytical func- 
tions (e.g. Sersic, exponential, Ferrer, Moffat, Gaussian) in 
multiple combinations. 

Principle in the available GALFIT fitting functions used 



throughout this paper is the Sersic profile (Sersic 1963 1968 



[Graham &: Driver|2005| which describes how the galaxy light 
profile varies as a function of radius. The Sersic equation 
provides the intensity 7 at a given radius r as given by: 



/ (r) = h exp 



l/n 



(1) 



where Je is the intensity at the effective radius , the radius 
containing half of the total light, and n is the Sersic index 
which determines the shape of the light profile (see Figure 
[2J. The value of &„ is a function of Sersic index and is such 
that r(2n) = 27(2n, h,,f\ where F and 7 represent the com- 
plete and incomplete gamma functions respectively! Ciotti 
1991 1. Varying the Sersic index parameter n allows one to 



model a wide range of galaxy profile shapes, with n = 0.5 
giving a Gaussian profile, n = 1 an exponential profile suit- 
able for galactic disks, and n = 4 a de Vaucouleurs profile 
commonly associated with massive spheroidal components 
such as elliptical galaxies. 

The only inputs required by SIGMA are the imaging 
data itself and the locations of the primary galaxies within 
them which are to be modelled. All additional parameters 
and starting values for extra neighbouring objects in the 
field of view (secondary objects), including PSF evaluation, 
are determined by SIGMA on the fiy on a per-galaxy basis. 
All scripting and additional programming is written in the 
open source and freely available R programming language 
( R Development Core Team||2010 l. Further information on 
the invocation of SIGMA may be found in Appendix [X] 

SIGMA operates in a semi-modular fashion, with an 
overarching master script calling and linking several key 
modules within. Each module is specialised in performing 
and handling a different task. A summary of each module 
and its purpose is shown in Table |2j and a schematic of the 
SIGMA data-fiow process is shown in Figure [3] The average 
run-time to profile a single primary object is 15 seconds per 
processoi[^sustained over several hundred thousand objects. 

^ bn can trivially be calculated within R using the relation 6„ = 
qgamma(0.5, 2n), where qgamma is the quantile function for the 
Gamma distribution. 

''' Using current computer hardware at the University of St An- 
drews. This consists of a 16 core Intel Xeon E5520 server with 
48GB RAM. 




r/r. 



Figure 2. The Sersic profile (Equation[T]l describes how a galaxy 
light profile varies as a function of radius, shown here for five dis- 
tinct values of Sersic index n. (top) Surface brightness at a given 
radius, (middle) Flux contained within a given radius, (bottom) 
The magnitude offset between the total magnitude of the galaxy 
and the Sersic magnitude at a given radius. 



3.1 SIGMA Master Script 



When initialising SIGMA a number of options are specified. 
In addition to several expected inputs such as band, nam- 
ing conventions and the combination and types of GALFIT 
functions to be modelled, SIGMA also allows the workload 
to be split over several processors. Using a greater number 
of processors directly decreases the amount of time required 
to analyse every primary object in the input catalogue, and 
is limited only by available hardware. A full summary of 
SIGMA input options is given in Appendix \K\ 

To begin, SIGMA's master script loads into memory the 
entirety of the GAMA input catalogue and defines a naming 
convention for each primary object based on its own unique 
identifier (SIGMAJNDEX). A template master CSV cata- 
logue is created into which all of the output data will accu- 
mulate as SIGMA loops across each primary galaxy. Once 
the setup is complete, the master script will loop across every 
primary object in turn. If for any reason a primary galaxy 
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Figure 3. A flowchart describing SIGMA's operation. Required inputs are imaginf 
maps), and an input catalogue with a list of galaxy coordinates. 



data in the form of GAMA mosaics (including weight 



causes a software crash, with attempted fixes as detailed in 
subsequent sections unsuccessful, SIGMA will report how 
far it was able to progress and record a NULL result before 
proceeding on to the next primary galaxy in the input cata- 
logue. We now discuss each module from Table [2] in turn. 



3.2 cutterpipe: Image Cutout and Preparation 

The CUTTERPIPE module creates and prepares the fitting im- 
age to be fed into GALFIT. Version 2 mosaics of the three 
GAMA regions are used as an input to cutterpipe, with 
a full description of the construction and manipulation of 



these files given in Hill et al. (2011 1 and summarised in Sec- 
tion [2] 

CUTTERPIPe's first task is to create the core cutout of 
the science image and its associated weight map. Using the 



WCS information stored in the FITS header of the appropri- 
ate mosaiced image, cutterpipe converts the input RA/- 
DEC into an x/y pixel coordinate. The upper and lower 
limits of a 1201 x 1201 pixel (~ 400" x 400") region centred 
on the primary galaxy are determined. Using the NASA 
HEASARC package's CFITSIO subroutine library, namely 
the routine fitscopy, cutouts centred on the primary galaxy 
on both the mosaiced science image, swpim, and mosaiced 
weight map, swpwt, are created. These cutouts are named 
cutim and cutwt respectively, fitscopy was found to be the 
most efficient routine at dealing with the large mosaic files 
in use, able to quickly analyse the input file and read into 
memory only the relevant area of interest, thereby reducing 
file handling time significantly. 

The process of creating the GAMA mosaics alters a 
number of keywords in the FITS header in order to bet- 
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Module 


Description 


File 

Outputs 


CUTTERPIPE 
Section [0| 


Creates a science and weight map 
cutout from the master GAMA 
mosaics using the CFITSIO 
routine FITSCOPY (Figure and 
performs an additional local 
background sky subtraction using 


cutim 
cutwt 


STARPIPE 

Section |3.3| 


Determines which frames 
contributed flux to the primary 
galaxy, and creates a catalogue of 
stars that lie within these frames 
using Source Extractor. 


psfws 
psfwt 
psfct 


PSFPIPE 

Section [3^ 


Generates an empirical 2D PSF 
at the primary object position 
using PSF Extractor (Figure jojl. 


psfss 
psfim 
psfsr 


OBJECTPIPE 
Section |3.5| 


Calculates starting parameters for 
size, brightness, position angle 
and ellipticity for the primary 
galaxy and any secondary 
neighbours (galaxies and stars) 
using Source Extractor. A 
dynamic search algorithm is used 
to attempt to detect the primary 
galaxy. Elongated objects (such as 
satellite trails) are removed from 
the secondary catalogue and 
instead added to a bad pixel 


objct 
segim 


GALFITPIPE 
Section |3.6| 


Fits an analytical function in 2D 
to the science image using 
GALFIT. Both primary and 
secondary objects are modelled, 
with any detected errors/crashes 
flagged and a fix attempted on a 
per-galaxy basis (Figures [7] and 

i 


segfr 
extct 
objim 



rtot = 2 X rgg 



(2) 



Table 2. Summary of the modules that comprise SIGMA, a brief 
description of their purpose, and a list of the file outputs produced 
by each. 



ter describe the nature of the mosaiced data. The mosaic 
headers are copied over to cutim and cutwt during their 
creation. Several of these keywords are required later in the 
fitting process by GALFIT in order to generate a sigma-map 
(an image showing the la confidence interval at every pixel). 
CUTTERPIPE reverts these to typical pre-mosaic values which 
are more appropriate for a smaller single image rather than 
a larger mosaic. GAIN, RDNOISE, NCOMBINE and EXP- 
TIME are set to values of 0.5, 3, 1 and 1 respectivelj|^ 

An estimate of the local background sky is then made 
with Source Extractor (v2.8.6; [Bertin fc Arno uts|1996[) using 
a variable background grid in a 3 x 3 mesh configuration. 
Possible grid sizes are 32 x 32, 64 x 64 and 128 x 128 pixels. 
The size of the chosen background grid is dependent upon 
the size of the primary galaxy: larger galaxies will lead to a 
larger background grid being used so as not to contaminate 
the sky estimate with galaxy flux. An initial basic estimate 
of the total size of the primary is given by: 



° These typical values are averages taken from pre-mosaiced data 
frames. 



where rgg is the radius of the primary galaxy which contains 
99% of the flux. This is obtained from the Source Extractor 
parameter FLUX_RADIUS setting PHOT_FLUXFRAC= 
0.9!|^ If rtot < 128, cutterpipe rounds up vtot to the 
nearest available grid size and performs a background sub- 
traction on the science image as appropriate. If rtot 5S 
128, no background subtraction is necessary, as the mas- 
ter GAMA mosaics have a 256 x 256 pixel background grid 
subtraction already applied. Although the value of actual 
subtracted sky varies with position on the cutout image, 
the specific value at the position of the primary galaxy, psky 
is recorded through the Source Extractor parameter BACK- 
GROUND. The error on the background sky estimate is then 
given by: 



Csky 



(3) 



y/0.9 X Ux X Uy 

where asky is the RMS of background sky counts across the 
cutout, and tIj, and Uy are the dimensions of the cutout 
in the x and y dimensions respectively. The background sky 
typically encompasses ~ 90% of any given cutout, and hence 
a factor of 0.9 is introduced into the above calculation to 
account for this. After extensive testing, this variable back- 
ground mesh method for sky estimation was found to be 
the most robust at removing small-scale sky fluctuations in 
the data without subtracting real galactic light from objects 
(see Section 



4.2.1 



for sky results discussion). 
An example image cutout, weight map and background 
estimation map are shown in Figure [4] 



3.3 starpipe: Star Detection 

STARPIPE uses Source Extractor to create a catalogue of star- 
like objects with which to create a point-spread function 



(PSF) in the subsequent psfpipe module (section 3.4 1. 

The first step is to determine which of the original 
pre-mosaic frames contain the primary galaxy. This step 
is non-trivial, as a single cut-out image (cutim) may con- 
tain data from several pre-mosaiced frames overlapping at 
random angles to each other, with only some of the frames 
contributing fiux (and therefore seeing information) to the 
primary galaxy. Calculating frame ownership is crucial in 
PSF determination, as using stars from non-contributing 
frames would skew the PSF estimate away from its true 
shape at the position of the primary galaxy. A method was 
devised to determine contributing frames using the inform- 
ation within the GAMA weight maps. Each pre-mosaiced 
frame is assigned a numerical value based on the global 
variance of the data for that frame. This value, repeated 
for each pixel, becomes the weight-map for that individual 
frame. Weight-map values are essentially unique to several 



^ Elsewhere in this paper, Source Extractor FLUX RADIUS will 

typically refer to re, a radius containing 50% of the flux of the 
primary galaxy. It is worth noting however that a size estimate 
produced by Source Extractor in this manner is known to be 
smaller than the true galaxy value, scaling as a function of Sersic 
index and thus absolute magnitude. This effect has been accoun- 
ted for, and does not adversely affect any of the analysis or results 
presented in this paper. 
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Figure 4. A series of setup images centred on G00196053. a) A science image cutout of 1201 X 1201 pixels (^^ 400" X 400") in the r band, 
with the primary galaxy centred at 601 X 601. The image is scaled logarithmically from —1 cTsky to 40 crskyj where crgky is the typical RMS 
of the sky in the r band. 6) A weight map cutout of 1201 X 1201 pixels (~ 400" X 400") in the r band, with the primary galaxy located 
in an overlap region at the white cross-hairs. Darker shades indicate larger values, therefore a greater number of overlapping frames, c) 
The additional local background sky removed from this cutout. The z-axis range shown is ±0.5 counts. For this galaxy, a background 
mesh size of 128 X 128 pixels (white square) was chosen in addition to the coarse 256 X 256 background mesh already applied in creating 
the GAMA mosaics. 



significant figures, and therefore useful in identifying that 
particular frame. During the SWarp process, overlapping 
imaging data is median combined (setting the SWarp ar- 
gument COMBINE_TYPE to MEDIAN) whereas weight- 
maps are co-added to produce a global weight-map repres- 
enting the change in the variance across the data. When 
two or more frames overlap, their individual weight map 
values are summed. Larger values indicate a greater number 
of overlapping frames. 

The value of the weight map at the primary position is 
determined, with all pixels of that value clearly contributing 
flux to the primary galaxy. This defines the initial primary 
region. However, since this primary region may be an overlap 
region itself, parent frames must also be determined. The 
weight-map values of all bordering pixels to the primary 
region are determined. Higher pixel values indicate a region 
which contains data from additional frames that did not 
contribute flux to the primary region, and so these pixels 
are discarded. Lower values (if any) indicate parent frames of 
the primary region, and (if they exist), their pixel positions 
are added to the primary region. This process will continue 
until all pixels are accounted for across the cutout. As an 
example of frame determination, contrast Figure[4]3 with the 
shaded red regions in Figure [7] 

Pixel determination via this technique is time intensive 
for the full 1201 x 1201 cutout region, as it requires analysis 
of 1.4 million pixels for each galaxy. A more efficient method 
is to reduce the number of pixels that require analysis by 
simplifying the weight map to its minimal number of pixels 
which still describe the nature of the data. Duplicate rows 
and columns in the weight map are removed, producing a 
simplified weight map, psfws, typically of order ~ 100 x 100 
pixels. This thereby reduces the number of pixels needing 
to be analysed by a factor of ~ 150, significantly speeding 
up the primary region determination. Figure |5] is one such 
simplification of the cutout weight map shown in Figure |4]3, 
in this case reducing the number of pixels from 1.4 miUion 
to ~ 20, 000. 




"I 1 1 1 1 1 r 

20 40 60 80 100 120 140 

Figure 5. A simplified version of the weight map shown in Fig- 
ure |4|). A single cutout from a GAMA mosaic may contain data 
from several frames observed on different nights in different seeing 
conditions. The weight map allows us to determine which frames 
contributed flux to the primary galaxy and which did not, thereby 
allowing a PSF to be constructed representing the seeing for that 
galaxy. Simplification of the weight map removes superfluous in- 
formation contained within the original cutout weight map, re- 
ducing the number of pixels to analyse typically by a factor of 
~ 150, and significantly speeding up starpipe. 



Once a primary region is determined, a local star cata- 
logue must be created. A modified version of cutwt, psfwt, is 
created, setting all non-primary pixels to a weight value of 
zero. This will bar Source Extractor detecting any objects 
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in those regions. The central 25 x 25 pixel region is also 
set to a weight value of zero to ensure against the primary 
galaxy being falsely classified as a star and used in later PSF 
analysis. 

A Source Extractor parameter file is created to out- 
put NUMBER, XJMAGE, YJMAGE, FLUX_RADIUS, 
FLAGS, FLUX_APER(1), FLUX_MAX, ELONGATION, 
VIGNET(25,25) and BACKGROUND. The numeric values 
following VIGNET determine the ultimate size in pixels of 
the 2D PSF created subsequently in psfpipe. A detection 
threshold of 2a above the background is specified, along with 
a minimum object detection area of 10 pixels, a 25, 000 count 
saturation level and a fixed sky pedestal value of zero counts. 
The image is filtered through a 5 x 5 pixel Gaussian convolu- 
tion kernel with F = 2 pixels. Source Extractor defaults are 
used everywhere else. Using these settings. Source Extractor 
is run across the cutout image cutim using the weight map 
image psfwt. An output catalogue, psfct, is created in the 
FITS_LDAC format, a format which saves the image data 
as well as the catalogue. 

3.4 psfpipe: PSF Estimation and Creation 

The point spread function (PSF) describes the blurring ef- 
fect of both the atmosphere and the telescope optics on our 
imaging data. Observed galaxy images have had their fiux 
redistributed according to this PSF. The galaxy fiux most 
affected by the PSF blurring is that which emanates from 
the core regions, where the gradient of the light profile is 
at its steepest. It is therefore crucial to have a good under- 
standing of the PSF when considering fitting smooth ana- 
lytical galaxy models to imaging data. Furthermore, most 
current galaxy modelling software weights model fitting to- 
wards higher signal-to-noise regions (typically the same core 
regions), increasing the importance of accurate, reliable PSF 
estimation. 

The PSFPIPE module is a wrapper around the PSF 
extraction software PSF Extractor (PSFEx v3.3.4; Bertin, 
priv. commp^, and produces a 2D PSF model to be taken 
into account at the later galaxy modelling stage. PSFEx ex- 
tracts precise models of the PSF from images pre-processed 
by Source Extractor, allowing for a wide range of PSF's to 
be quickly and accurately constructed, including arbitrary 
non-parametric features present in the PSF. 

In brief, the sample of objects from the psfct catalogue 
created by STARPIPE is initially used for analysis. PSFEx 
reduces this object list to a star sample based on a set of 
pre-defined criteria. A signal-to-noise limit of at least 10 is 

required, and objects with an eccentricity of ^ \a+b) ) ^ ^'^^ 
are removed, where a and b refer to the semi-major and 
semi-minor axes respectiveljj^ Each star's full-width-half- 
maximum, r, is estimated, with only stars in the pixel range 

More information on the PSFEx software may be found at 
http:/ /www. astromatic.net/software/psfex. 

PSFEx refers to this quantity as elhpticity rather than eccent- 
ricity, however its definition is more akin to that of the latter. 
We adopt the terminology eccentricity here to avoid confusion 
with the standard definition of ellipticity used throughout the re- 
mainder of this paper, namely e = 1 — -. An eccentricity of 0.05 
therefore corresponds to an ellipticity of e ~ 0.095. 



2 < F < 10 accepted. Furthermore, variability in the star 
sample is limited to the central 50% quantile. After extensive 
testing on the variation in PSF quality with star sample 
size, and communication with the authors of PSFEx, we 
found that a star sample size of at least 10 stars is necessary 
to ensure that the resultant PSF is not adversely affected 
by small-number biases. Therefore, if fewer than 10 stars 
remain in the star sample after selection criteria have been 
applied, SIGMA will loop back to CUTTERPIPE and expand 
the cutout region to 1501 x 1501 pixels (~ 500" x 500"). 
The mean number of stars used for PSF estimation in the 
r band is 24.4, with 10.2% of cutouts containing fewer than 
10 stars after the cutout region has been expanded. 

Cutout images of each star are pre-stored in the 
FITS_LDAC format of psfct, the size of the cutout having 
been specified at the Source Extraction stage. PSFEx uses 
the positional information from Source Extractor to mask 
nearby neighbours to the final star sample, and presents this 
sample in the output psfss FITS image (Figure |6|i) . 

The variation in the shape of the stars in the star sample 
is then modelled in both x and y as a function of position 
in the field by a 2D n**^ order polynomial function. Higher 
order terms in the fit (i.e. x, x^, etc.) describe the variation 
in the PSF at positions away from the centre of the frame. 
Within SIGMA, the primary galaxy is always centred in the 
cutout image, and so a zeroth order polynomial was found 
to adequately describe the PSF. The best-fit polynomial is 
sampled at a 1:1 ratio relative to the input data, and an 
output PSF image psfim is produced of the same size as the 
input cutout stars, 25 x 25 pixels (Figure [Hja). As a con- 
sistency check, scaled models of psfim are fit to each of the 
input stars in psfss, and a residual map psfsr produced (Fig- 
ure [6|;). Note that some of the PSF residuals still show no- 
ticeable structure once the PSF model has been subtracted 
from the star sample. This is as expected when subtracting 
a zeroth order PSF model (only accurate at the location of 
the primary galaxy in the centre) from a star sample taken 
over a large area on the sky. Those stars with noticeable 
residuals therefore are typically either significantly spatially 
separated in the field of view from the primary galaxy or 
approaching saturation (or both) . Both of these factors are 
accounted for by PSFEx when constructing the model PSF. 
The stars chosen as part of the star sample are shown in 
orange circles in Figure [T] with each circle numbered ac- 
cording to their position in Figure [6^, starting at 1 in the 
bottom left and increasing horizontally left-to-right and then 
bottom-to-top. 

3.5 objectpipe: Object Detection 

A second catalogue of objects optimised for galaxy detec- 
tion, objct, is created in objectpipe, to be later fed into 
GALFITPIPE. This catalogue provides the basic starting 
parameters for the primary galaxy and any secondary 
galaxies and stars in the frame, objectpipe also creates a 
segmentation map of the frame, modified to mask any erro- 
neous regions of flux in the image which may cause fitting 
problems (e.g., satellite trails). A Source Extractor para- 
meter file is created containing X_IMAGE, Y_IMAGE, 
MAG_AUTO, FLUX_RADIUS, KRON_RADIUS, 
AJMAGE, BJMAGE, THETAJMAGE, ELLIPTICITY 
and CLASS STAR. These give position (xjy), luminosity. 
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Figure 6. PSFEx generates an empirical PSF (b) from a host sample of representative stars (a). This figure shows 63 sample star 
cutouts of 25 X 25 pixels each chosen from around GAMA object G00196053, whose real sky positions may be noted in Figure [t] (orange 
circles). Panel (c) represents the residual of each star sample with a scaled form of the PSF subtracted from each, (a) and (c) are scaled 
logarithmically from —1 Ugky to 40 (Tskyi where (Tsky is the typical RMS of the sky in the r band. 



size, position angle and ellipticity for the primary and 
all secondaries in the field. Source Extractor settings 
are similar to those used in STARPIPE, excepting a lower 
detection threshold of 1.8cr above the background (where 
0" is the RMS as estimated by Source Extractor), and a 
lower deblending contrast parameter of 0.0001. The Source 
Extractor Neural-Network Weights VI. 3 file is used in 

the creation of the CLASS STAR parameter, as well as 

a standard 5x5 pixel Gaussian filter, F = 2 pixels, used 
during object detection. 

OBJECTPIPE calls Source Extractor, and records the res- 
ults. If initially the primary galaxy is unable to be located 
within a 5 pixel radius of the input coordinates, OBJECTPIPE 
will in the first instance decrease the detection threshold in 
steps of 0.4a down to Icr above the background until an 
object is found, re-running Source Extractor as appropri- 
ate. This usually occurs with faint objects in the field, or in 
crowded regions. If the primary object is still unable to be 
located, the threshold is reset to 1.8a, and a larger search ra- 
dius of up to 15 pixels from the input coordinates in 5 pixel 
steps is tried. This stage accounts for large nearby galaxies 
whose centroids are not matched to better than 5 pixels, 
hence requiring a larger detection area. If multiple-matches 
are found, the largest object will be taken to be the primary 
galaxy. If at this stage the primary galaxy is still not found, 
OBJECTPIPE will report a null detection, and move on to the 
next primary in the input catalogue. 

Output parameters from Source Extractor are modified 
by OBJECTPIPE before being fed into GALFITPIPE with the 
exception of magnitude which is used unaltered. Position 
angle is modified to the GALFIT standard (by adding 90 
degrees), increasing anti-clockwise from the positive x axis. 
Ellipticity e is converted to an axis ratio using the relation: 



(0.32 X r2 



(5) 



e = 1 - - 
a 



(4) 



where rso is the (unmodified) Source Extractor half- 
light radius as given by FLUX_RADIUS (setting 
PHOT_FLUXFRAC= 0.5) and F is the Full-Width Half- 
Maximum of the PSF of the primary galaxy. A minimum 
bound on re of 1 pixel is enforced. This conversion cor- 
rects for the fact that Source Extractor's output half-light 
radii are circularised and based on PSF convolved data, 
whereas GALFIT radii are along the semi-major axis and 
intrinsic (non-PSF convolved). The value of 0.32 was derived 
from simulated test data, see Appendix A of [Driver et al.| 
(2005) for further details. Figure [bT] shows a comparison of 
corrected and uncorrected radii against modelled GALFIT 
radii for all GAMA objects in the r band. This suggests 
that the revised starting value for is appropriate and an 
accurate first estimate of the true half-light radius of the 
primary galaxy. Due to the downhill minimisation employed 
by GALFIT, it is important to provide input parameters as 
close as possible to the desired solution in order to avoid 
local minima. 

Once physical parameters for the primary galaxy have 
been determined, a segmentation map of the frame is cre- 
ated to be used as a potential mask for secondary features 
should modelling them fail. Secondary objects whose el- 
lipticity is greater than 0.95 are excluded from modelling 
and will instead be masked, as these are determined to be 
satellite trails or bad data, and consequently difficult to 
model. Similarly, secondary objects with a stellaricity index 



with semi-minor axis h and semi-major axis a. Half-light 
radius re is estimated using the relation: 



of CLASS_STAR> 0.8 (see [Bertin fc Arnouts|[l996l l are 
modelled by a PSF within GALFITPIPE, with all others be- 
ing modelled using a single Sersic function. A relatively low 

CLASS STAR boundary is chosen as tests have shown that 

a more reasonable fit is produced when fitting a PSF rather 
than a Sersic function to ambiguous objects. A graphical 
representation of detected galaxies, stars, weight-map areas 
and secondary neighbour determination for GOOl 96053 is 
shown in Figure [7] 
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and mtot- For additional information on the operation of 
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GALFIT, refer to Peng et al. (2010i 
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Figure 7. Final detail analysis plot for G00196053. Green shaded 
areas represent all detected objects, with stars chosen as part of 
the PSF star sample circled in orange. The GALFIT fitting area 
is outlined in yellow, and the weight map frame edges in blue. 
Secondary objects which will be modelled as nearby neighbours 
in the fitting process are represented in turquoise. Note how no 
sample stars are taken from frames which did not contribute dir- 
ectly to the flux of the primary galaxy (areas represented with 
red shading). Consequently, no stars are taken for PSF analysis 
from the red shaded areas. The image is scaled logarithmically 
from —1 (Tsky to 40 (Tgijy, where a-gity is the typical RMS of the 
sky in the r band. 



3.6 galfitpipe: Galaxy Fitting 

The actual modelling is handled by galfitpipe, which is 
a wrapper around the GALFIT image analysis software 
(v3.0.2; Peng et al.|2010 1 along with several event handlers 
and logical filters written in R. GALFIT is a 2D parametric 
galaxy fitting algorithm written in the C language. It allows 
for multiple parametric functions (such as Sersic, exponen- 
tial, Ferrer, Moffat, Gaussian, etc.) to be modelled simul- 
taneously as either multiple components of a single object, 
multiple objects in a single frame, or combinations thereof. 

GALFIT uses a Levenberg-Marquardt algorithm to fit 
a 2D function to 2D data, in doing so minimising the global 
until the gradient A^^ has become negligible and conver- 
gence is reached. When a global minimum is thought to be 
found, GALFIT introduces a 10 iteration cool-down period, 
sampling the parameter space around the best-fit paramet- 
ers in an attempt to overcome the problem of converging on 
a local rather than global minimum. 

In this paper we fit each primary galaxy with a single 
Sersic function containing 7 free parameters: object centres 
xo and j/o; total integrated magnitude mtot', effective radius 
along the semi-major axis re; Sersic index n; ellipticity e 
and position angle 6. Secondaries (galaxies and stars) will 
also be modelled by either a Sersic function or a scaled PSF 
as appropriate. The PSF contains 3 free parameters; xq, yo 



All primary inputs to GALFIT are taken from Source 
Extractor and modified as described in Section |3.3[ with 
the exception of Sersic index which starts at ninitiai = 2.5. 
After extensive testing on the r band data, it was found 
that the chosen starting Sersic index has little to no effect 
on the end result, and so choosing a value in the centre 
of the expected parameter range was deemed appropriate. 
No explicit constraints were put on the range of acceptable 
Sersic indices upon which GALFIT may converge, however, 
GALFIT has internal limits of 0.05 < n < 20, where the 
lower limit is a 'soft' limit (indices scatter around this value) 
and the upper is a 'hard' limit (indices may not converge 
above this value) . More conservative limits were not enforced 
on Sersic index so as not to lead the final results and make 
presumptions about Sersic index distributions. More detail 
on chosen initial conditions may be found in Appendix [B] 

In order for galfitpipe to function correctly, it needs 
the cutout science image from CUTTERPIPE, cutim; the asso- 
ciated segmentation map and object catalogue from object- 
pipe, segim and ohjct respectively; and a 2D FITS image 
PSF representing the PSF at the primary galaxy location 
from STARPIPE and psfpipe, psfim. Note that the weight 
map (cutwt) is no-longer required at this modelling stage. 

Once the aforementioned files are in place, an initial 
fitting region radius on the cutout is defined by: 



= 2rKro„ (1 cos {9) ] + (1 - e) | sin {6) \) 
Ty = 2r-K.on (I sin (0) \ + {l-e)\ cos (0) |) 



(6) 
(7) 



in order to account for the ellipticity e of the object and its 
position angle 0. Objects within the central 2rx x 2ry of the 
fitting region will be convolved with the supplied PSF at 
the modelling stage. The segmentation map is modified to 
unmask all secondary objects in the fitting region, with the 
resultant map saved into a new segfr file. 

A GALFIT feedme file is created containing the starting 
values for every object being modelled (primary and second- 



ary) as described in Section 3.3 and above. A constraints 



file is used to constrain secondary objects. These objects 
are constrained in order to reduce the fitting time, and re- 
duce the size of the allowed parameter space, xq and yo are 
constrained to ±3 pixels of their input parameters, ellipti- 
city is constrained to < e < 0.95 and half-light radius is 
constrained to re,initiai/4 < re, final < 4re, initial. A final para- 
meter for sky is added to the bottom of the GALFIT feedme 
file, fixing the value of the sky to zero counts. 

GALFIT is then initialised, fixing the sky RMS to that 
measured in OBJECTPIPE. The time taken to converge on a 
fit scales with the size of the fitting region, and the number 
of secondaries being fit. Once the GALFIT process has fin- 
ished, its output (if any) is read and processed, galfitpipe 
scans the primary galaxy for a number of problems in this 
order: 

(i) Crash or a segmentation fault 

(ii) Galaxy centre migration of -y/x^ + j/^ > r^ initial 

(iii) An exceptionally large radius of logio ^ ^ > 3 

(iv) An exceptionally small radius of logio ^ J^";^""'^ ^ < 3 

(v) A high ellipticity of e > 0.95 
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If any of these are detected, a fix will be attempted and 
GALFIT re-run as appropriate. Fixes attempted vary de- 
pending on the problem encountered. If a crash or segment- 
ation fault are detected, GALFIT will be re-run modelling 
only the primary galaxy, with all secondaries masked. This 
usually occurs for large nearby objects with a high number 
of secondary neighbours and foreground stars, providing the 
Levenberg-Marquart minimisation routine in GALFIT with 
many local-minima. If the centre migrates away to fit a sec- 
ondary feature, GALFIT will be re-run with the primary 
centroids fixed to their starting values. Large or small radii 
are initially handled by suggesting a lower starting shape 
parameter (n = 0.5). This usually assists GALFIT in finding 
a way out of any local minima. If this attempt still provides 
a wildly different size to the input parameter, the size is 
fixed to the input and GALFIT re-run. Finally, a high el- 
lipticity usually indicates the model has migrated away to 
fit flocculent secondary features. Re-running GALFIT with 
a starting ellipticity of e = 0.1, i.e., highly circular, in most 
cases mitigates this problem. If all fixes have been attemp- 
ted and problems persist, GALFITPIPE will record GALFIT's 
best-guess model parameters and move on to the next ob- 
ject in its catalogue, with a flag updated to reflect the fitting 
history. 

If the fit has been successful, the output multi-HDU 
FITS file from GALFIT is saved as objim, with a catalogue 
of final modelled secondary objects saved as extct. An ex- 
ample model output for GAMA galaxy G00092907 is shown 
in Figure [8] A series of value added measurements are cal- 
culated and added to the structural measurements already 
taken for the primary galaxy. These include /io (central sur- 
face brightness), (surface brightness at the half-light ra- 
dius), (fie) (average surface brightness within the half-light 
radius) and rgo (radius along the semi-major axis that con- 
tains 90% of the tota|^ Sersic flux), amongst others. These 
values, along with the output parameters from GALFIT and 
previous SIGMA modules are added to a comma-separated 
variable (CSV) catalogue, allowing SIGMA to move to the 
next primary galaxy in the input catalogue. 



4 SIGMA OUTPUT 

The SIGMA master catalogue, SersicCatv?, provides meas- 
urements of Sersic index, half-light radius, position angle, 
ellipticity and magnitude in addition to extra pre-modelling 
sky estimation. Source Extraction and PSF Extraction 
measurements and post-modelling value added measure- 
ments as detailed in Section [3] Magnitudes contained within 
this catalogue are deflned according to the AB magnitude 
system, and have not been corrected for the effects of fore- 
ground Milky Way dust extinction. The catalogue is an out- 
put of the GAMA SersicPhotometry data management unit, 
and contains 527 columns of data; 58 columns per passband, 
and 5 additional common descriptive columns. Here we dis- 
cuss the results of the modelling pipeline. 



Name 


Number 


Definition 


TilingCatvll 


169, 850 


Complete GAMA tiling 
catalogue 


ScrsicCatvOr 


167, 600 


Removes star-like objects 


Survey 


150,633 


Removes rpetro < 19-8 


Common 


138,269 


Requires SIGMA Source 
Extractor coverage in 
{ugriz)+Y+J+H+K 


Matciied 


116,951 


Requires a match in the 
StellarMassesv03 catalogue 



Table 3. Table defining various sample definitions in use through- 
out this paper. Cuts are sequential, and include the definitions 
from previous rows. 



4.1 Sample Definitions 

Table [3] summarises the various sample definitions in use 
throughout this paper. Our initial input is the GAMA tiling 
catalogue, TilingCatvll, which contains 169, 850 sources. Of 
these sources, 167, 600 are classed as galaxy-like (as defined 
Baldry et aL]|20Tol l. SIGMA was run across this galaxy- 



^"^ integrated to infinity 
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like subset independently in all nine bands, the output of 
which defines the SIGMA master catalogue, SersicCatvl, 
which is available via the GAMA database. However, we 
define additional sub-samples in order to facilitate further 
analysis of the data throughout the remainder of this paper. 
This is to ensure that selection bias does not adversely af- 
fect our conclusions. SersicCatvOl contains sources fainter 
than the deepest nominal GAMA limit of rpetro = 19.8, 
and so a cut was made limiting sources to rpetro < 19.8. 
A common coverage sample was constructed so as not to 
compare galaxies between bands whose observations are 
incomplete or have missing data. This was defined using 

the Source Extractor Auto magnitude SEX MAG X from 

SersicCatvl, a product of the OBJECTPIPE module, where 
X=UGRIZYJHK. A common region is defined as hav- 
ing a detected Source Extractor magnitude in any of the 
SDSS bands (ugriz) as well as in each of the UKIDSS bands 
{YJHK). Incompleteness mainly affects the NIR bands, 
with noticeable UKIDSS footprint gaps visible in the final 
common coverage area shown in Figure [l] The number of 
detected sources in individual SDSS bands is typically very 
high, > 97%, with the exception of the u band. The u band 
data has a detection percentage via this method of 50.8%, 
indicating the poorer quality of the data in that band. For 
this reason, u band data is excluded from further fits to the 
data, with relations instead extrapolating into the u band 
wavelength for reference. The common coverage area reduces 
the sample to 138, 269 galaxies and is used throughout Sec- 
tion |4]and the beginning of Section [5] (with the exception of 
Figure [l7| . A full listing of detected and modelled sources 
used in each band may be found in Table |4] In any analysis 
that makes use of data from the StellarMassesvOS catalogue 
(rest-frame K-corrected u — r colours or stellar masses), 
a reduced matched coverage sample of 116,951 objects is 
defined. 
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Figure 8. An example model output for GAMA galaxy G00092907. The original SDSS r band image is shown in the top left. The 
2D model of this galaxy and its residual (image - model) are also shown as indicated. Inset into the residual is a postage stamp of 
the PSF constructed for this galaxy. Blue captions within the PSF postage stamp indicate (anti-clockwise from top-left) the full-width 

2 

half-maximum of the PSF; the number of stars used in creating the PSF; the ^ for the PSF model fit and the size of the PSF postage 
stamp (8.5" X 8.5"). The ID profile for this galaxy, calculated by taking the average counts along ellipses centred on the primary galaxy 
and displayed against the semi-major axis, was created using the IRAF package ELLIPSE. Example ellipses from ELLIPSE have been added 
to the image in the bottom left to guide the eye, spaced evenly at intervals of 2" along the semi-major axis. Note that . Inset into the 

2 

ID profile are relevant output Sersic modelling parameters for this galaxy, and the overall ^ for the Sersic model is shown in grey at 
the bottom of the figure. The ID profile is a ID measure of 2D data. Any flux from secondary objects (neighbouring galaxies and stars) 
lying outside the ID mask and overlapping with the primary galaxy will be counted as belonging to the primary galaxy by ELLIPSE. The 
ID profile therefore should chiefly be used as a guide as to the true light distributions of both the 2D image and 2D model. For plotting, 
the image data is divided by some scaling constant (150 counts in the r band) and scaled using the arctan function with cut levels at 
-f and f . 





Detected 


Modelled 


u 


85, 138 


81,120 


9 


165,367 


165,196 


r 


166,506 


166,384 


i 


166,675 


166,377 


z 


163,902 


160,684 


Y 


156,702 


156,280 


J 


152,316 


151,612 


H 


159,464 


158,797 


K 


157,537 


156,662 



Table 4. The number of detected and modelled galaxies in Ser- 
sicCatv07 for each band. SersicCatvOl contains 167,600 galaxies 
in total. 

4.2 Analysis 

4-2.1 Additional Sky Subtraction 

As part of the cutout creation process, SIGMA uses a vari- 
able background mesh to estimate and subtract the back- 
ground sky for each galaxy in each band before any other 
image analysis takes place. Sky correction distributions are 
mostly Gaussian in shape, with a small bias to recovering 
positive sky values most likely owing to background source 
contamination at the sky estimation stage. The additional 



correction on top of that already applied at the GAMA mo- 
saicing stage is usually small. In the r band for example, 
the sky correction distribution has a 3(T-clipped mean of 
0.56 ADU's (~ 0.01% of the sky pedestal) and a standard 
deviation of cr = 2.06 ADU's. Longer wavelengths produce 
larger corrections as expected. The accuracy with which we 
were able to estimate the sky using this method was found 
to produce good quality sky estimates in an efficient and 
relatively fast manner. 

An additional spike feature at zero counts relates to ob- 
jects whose determined preferred background mesh size was 
larger than that already used in the creation of the GAMA 
mosaics. If this has occurred, SIGMA performs no sky sub- 
traction, and returns zero counts. This feature affects 0.39% 
of galaxies in the r band, and so whilst a larger mosaicing 
background mesh may be preferred for future surveys, it is 
not believed to be a major issue affecting these data. 

4-2.2 Astrometry 

Initial checks were made on the output astrometric accuracy 
of the SIGMA models. Figure [9] shows the computed astro- 
metric offsets between the input SDSS positions and their 
GALFIT modelled SIGMA positions for all bands, with each 
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Figure 9. Astrometric offsets in RA and Dec between the input SDSS positions and the modelled SIGMA positions for all nine bands. 
Contours range from the 10*'' to the 90*'' percentile in steps of 10% with the peak density of each distribution represented by a yellow 
cross. Each sub-plot is exactly 2x2 pixels in dimension. The global systematic offsets in the NIR UKIDSS data (YJHK), typically 
~ 0.07" (~ 0.2 pixels), are caused by minor variations in the WCS definitions between SDSS and UKIDSS. SIGMA accounts for this 
during the modelling phase. 



sub-plot representing 2x2 pixels. Generally speaking, the as- 
trometry is in good agreement with SDSS, with the r band 
offset Cr = 0.010" (0.029 pixels) and a 1-sigma spread of 
lar = 0.044" (0.130 pixels). The one exception to this is the 
u band data, showing a much larger spread in the recovered 
centroids owing to the poor quality and depth of the data in 
this band. There are however two interesting features worthy 
of note in this figure. First, the apparent asymmetry in the 
SDSS astrometry, particularly noticeable in the higher qual- 



ity r and i bands. Second, the global systematic offsets in 
the NIR UKIDSS bands (YJHK) of approximately 0.07". 

The asymmetry present in the SDSS astrometric data 
is found to be associated with an individual SDSS strin^'^ 



In SDSS imaging data, a single run covers a strip. Two strips 
constitute a stripe, with the second strip offset from the first in 
order to cover a continuous area. 
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that crosses the G09 field at a large angle of incidence with 
respect to the equatorial plane. Galaxies whose input ima- 
ging data lie in this strip appear to have centroids scattered 
around ARA ~ -0.05", ADec ~ 0.10" rather than the 
origin. This feature is less prominent in the lower quality 
SDSS bands as it becomes lost in the random scatter, and 
consequently the effect is most noticeable in the r and i 
SDSS bands. Since this asymmetry affects only one strip 
of an SDSS stripe, the error must have been introduced at 
the splicing stage within the SDSS pipeline. These offsets 
remain small however, and are not believed to seriously af- 
fect this study as they are accounted for during the SIGMA 
modelling pipeline. 

Global systematic offsets in the NIR bands represent 
minor differences in the WCS calibration between the SDSS 
and UKIDSS data. Any discrepancy between the imaging 
data would be carried through to the larger GAMA mosa- 
ics. This feature also varies according to GAMA region, with 
measured offsets of approximately 0.05", 0.11" and 0.09" in 
G09, G12 and G15 respectively. As with the previous fea- 
ture, whilst consistent, offsets remain small sub-pixel vari- 
ations (~ 0.2 pixels) and therefore are not believed to be a 
major factor affecting cross-band matching between sources 
within GAMA. These features do not arise at the GALFIT 
modelling stage, as similar plots comparing input SDSS pos- 
itions against pre-modelling Source Extractor centroids from 
SIGMA show similar results albeit with larger spreads. On 
the contrary, SIGMA should do a better job of recover- 
ing true centroids due to GALFIT's model extrapolation 
method in estimating centroids. This makes SIGMA robust 
against astrometric errors such as this by recentring every 
galaxy at the modelling stage, emphasising the strength of 
full modelling against basic source extraction. 

4.2.3 Seeing 

An independent measure of the seeing and the form of the 
PSF for each galaxy in each band is a necessary require- 
ment when considering galaxy modelling. Through PSFEx 
in PSFPIPE, SIGMA is able to provide robust measurements 
of the PSF for each galaxy as described in Section |3.4| 
prior to the GALFIT modelling stage. Figure [To] shows the 
recovered PSF full-width half-maxima F for every galaxy 
within the SIGMA common sample for all 9 bands. Each 
density curve has a main peak in the range 0.7" < F < 1.4", 
and an additional peak at F = 0.4", which shall be discussed 
later. 

We note that on average the NIR data is of better see- 
ing than the optical, with the former in the range 0.6" < 
F < 1.3" and peaking at around F = 0.9", and the latter in 
the range 0.8" < F < 1.7" with variable peaks. These ranges 
are in good agreement with UKIDSS {K band) and SDSS (r 
band) seeing targets of Fukidss,k < 1-2" and Fsoss.r < 1-5" 
respectively. Some of the NIR data displays a secondary 
peak, particularly in the K band, possibly due to the use 
of microstepping in the taking of some of the UKIDSS data. 
The worst quality seeing data is in the u band, exhibiting 
the largest width distribution, and the highest seeing data 
on average. This distribution of its mean across the GAMA 
regions is represented in Figure [T] with the data points col- 
oured according to the measured seeing at that location. 
This figure shows significant striping in the SDSS r band 
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Figure 10. Recovered full-width half-maximum PSF values from 
the SIGMA common coverage sample. 

due to the drift scan mode of collection, and a measure of 
consistency coupled with lower average values across each of 
the UKIDSS bands. This could cause significant problems 
for image analysis routines, with average seeing doubling on 
the scale of a few pixels. Modelling the PSF and using that 
model at the galaxy modelling stage, as in SIGMA, goes 
some way towards mitigating this issue. 

An additional peak at F = 0.4" represents those frames 
where no stars were detected in order to compute the PSF in 
that region, and so a generic value of F = 0.4" is returned. 
Note that for the majority of bands this problem is minimal, 
becoming most noticeable in the lower quality u and 2 band 
data. 

4.2.4 Surface Brightness Limits 

Consideration of the surface brightness limit beyond which 
data becomes unreliable at the la level is also important. An 
estimate of the surface brightness limit at any given position 
may be given by 

Miim = ZP-2.51ogjRMS (8) 

where ZP is the zero-point of the imaging data, and /rms is 
the root-mean-square of the background sky per square arc- 
second. Note however that this provides a worst-case scen- 
ario value to the surface brightness limit, with the real limit 
likely to be deeper on a per-galaxy basis dependent upon 
the number of pixels (n) used in constructing the 2D model 
at large radii from the core region, and scaling as ^/n. Fig- 
ure [TTJshows the global distributions in for the SIGMA 
common coverage sample across each bandpass, with the me- 
dian surface brightness limits in each band given inset into 
the figure. We note that the shorter wavelengths typically 
exhibit deeper limits, as expected, with a transition occur- 
ring to shallower limits beyond the i-z interface. Figure [12] 
shows the spatial variation of /inm across the GAMA fields. 
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Figure 11. Apparent surface brightness limits for all galaxies 
within the SIGMA common coverage sample, with median surface 
brightness values for each band inset. 



The deepest ^lim data is represented by blue data points, 
the shallowest by red. The centroid weighting mechanism 
employed by GALFIT should minimise the impact of a spa- 
tially varying /ium, and therefore should not heavily affect 
the output results from SIGMA. 



4.3 Results 

4.3.1 Case Study Examples 

We present example model fits for an individual galaxy 
across all nine bands and various galaxies separated in mag- 
nitude space in Figures [13] and [14] respectively. These figures 
each represent the input 2D science image, model and re- 
sidual along with a ID profile radiating outwards from the 
core region of the galaxies along the semi-major axis. The 
displayed input image is a postage-stamp sub-region of the 
background corrected cutout from the original GAMA mo- 
saiced image. The yellow box in each of the input image 
postage stamps represents the size of the fitting region as 
determined in GALFITPIPE, the dimensions of which specify 
the size of the output model FITS image. Recovered Sersic 
parameters are listed inset into the ID profile plot. 

Figure [13] shows the output SIGMA results for the el- 
liptical galaxy G00032237 across all nine bands. Each image 
is modelled independently in each band, leading to a variable 
fitting region size dependent upon local conditions including 
object density and the physical size of the primary galaxy. 
The residuals in each band show the high quality of the 
fit for this particular galaxy, bar some minor core disturb- 
ance in the higher quality r and i bands. These bands cover 
wavelengths that are more sensitive to dust attenuation, in 
this case potentially highlighting small quantities of dust in 
the centre of the galaxy possibly related to a recent minor- 
merger or some form of morphological disturbance. Dust 



has the effect of perturbing the light profile slightly away 
from that of a purely single-Sersic object. Interestingly how- 
ever, no evidence for dust lanes are evident in the lowest 
wavelength u band residual. This should not be surprising 
considering the lower quality data of the u band, hence these 
small perturbations would be lost in the noise of the im- 
age. Barring the u band data, and despite dust attenuation, 
the recovered Sersic indices remain relatively stable, ranging 
from n = 4.21 to n = 4.50 in g to K. Sersic index peaks in 
the J band (rimax = 4.73) and reaches a minimum in the r 
band (nmin = 3.82, excluding u band data). Modelled ellipt- 
icity e and position angle 9 are also in good agreement, with 
recovered Sersic magnitude evolving as expected across this 
wavelength range. Interestingly, the recovered half-light radii 
show a size-wavelength dependence, ranging from r-g = 3.74" 
to Te = 2.13" in g to K. 

Secondary objects whose object centres lie outside the 
fitting region but whose flux leaks into it are masked so as 
not to effect the model fit. One such object can be seen 
in the upper-right corner of the fitting region in the r band 
postage stamp in Figure[T3] The GALFITPIPE module creates 
a bad-pixel mask using the segmentation map provided by 
OBJECTPIPE. Should GALFIT reach an error whilst trying to 
converge on a model, a potential additional fix is to mask all 
secondary objects in the field of view, and re-run GALFIT. 
This dynamic masking, whilst not the first choice for pro- 
ducing a model, typically allows high-quality and consistent 
model data to be extracted from difficult regions where oth- 
erwise none would exist. 

Figure [14] shows nine example galaxies in the r band 
from the SIGMA common coverage sample, separated in 
magnitude space in approximately equal SDSS r band mag- 
nitude steps of Arrir = 0.5 from the faintest GAMA limit 
of rrir = 19.8. These galaxies span a wide range of mor- 
phologies and environments, exhibiting the large variance 
in the input data processed by SIGMA. In each case, the 
residual images show the quality of the fits are relatively 
high, more so for obvious single-component objects such 
as the huge elliptical galaxy G00506119 (m^ = 15.8) than 
multi-component objects such as the barred-spiral galaxy 
G00369161 (rrir = 16.8). In the case of the latter, despite 
a global single-Sersic fit to a multi-component object, the 
resultant model does a good job at describing global para- 
meters such as magnitude. Despite a relatively disturbed fit 
to the secondary neighbour of G00177815 (nir = 18.3), the 
quality of the primary galaxy model remains high. This out- 
lines the exclusive use of secondary objects in accounting for 
additional flux in the wings of primary galaxies. Note that 
whilst the quality of model fits reflected through the resid- 
uals appears to become better at fainter magnitudes, this 
effect is more likely an example of non-resolved components 
of a galaxy. 

4.3.2 Global Results 

Complete distributions for the SIGMA common coverage 
sample of 138, 269 galaxies are shown in Figure [Ts] Here we 
plot ID density distributions for recovered model Sersic in- 
dices, half-light radii and ellipticities in each band. Alongside 
these distributions are displayed the average model galaxies 
based on median values from the aforementioned paramet- 
ers. 
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Figure 12. Apparent surface brightness limits for all galaxies within the SIGMA common coverage sample as a function of their position 
on the sky in right-ascension and declination. The three GAMA regions are displayed, as indicated, with each band labelled along the 
right side of the figure. Bands are offset in declination in order to differentiate them from one another. The K band is situated at the 
correct GAMA coordinates. Surface brightnesses are shown as offsets from the median surface brightness for that band, the values of 
which are found in Figure [TT] Blue data points represent the deepest limits and red the shallowest. 
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Figure 13. Model fits for G00032237 across all nine bands. Eacli column represents (from loft to right) the original input image, the 
model fit to the input image, the residual image (input - model) and the ID surface brightness profile along the semi-major axis (averaged 
along the annulus). The fitting region (the region within which 2D modelling takes place) is represented by a yellow box. Recovered 2D 
intrinsic (i.e., prior to PSF convolution) Sersic parameters are listed inset into the ID profile plot. The images are scaled logarithmically 
from —1 iTsity to no"s]jy, where a^y^y is the typical RMS of the sky in that band, and n is some scaling constant (generally, n ~ 40). 
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Figure 14. Model fits for nine galaxies in the r band separated in magnitude space by approximately Am^ = 0.5 in the range 
15.8 < rrir < 19.8. Each column represents (from left to right) the original input image, the model fit to the input image, the residual 
image (input - model) and the ID surface brightness profile along the semi-major axis (averaged along the annulus). If the fitting region 
(the region within which 2D modelling takes place) lies within the image thumbnail above, it is represented by a yellow box. If no yellow 
box is visible, the fitting region is larger than the thumbnail. Recovered 2D intrinsic (i.e., prior to PSF convolution) Sersic parameters 
are listed inset into the ID profile plot. The images are scaled logarithmically from —1 o-gky to 40 iTsky> where asiiy is the typical RMS 
of the sky in the r band. 
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Recovered Sersic parameters peak primarily in the 
range 0.2 < n < 10, with additional peaks at n ~ 0.05 
and n = 20 arising due to failed fits, discussed in more de- 
tail below. The primary range appears bimodal in nature, 
consisting of two approximately Gaussian-like distributions 
whose means are n ~ 1 and n ~ 3.5. These two peaks, for 
the most massive/brightest systems in GAMA, correspond 
to the two main galaxy morphologies as originally identi- 
fied by Hubble, namely, late-type disk-dominated galaxies 
and early-type spheroid-dominated galaxies for n = 1 and 
n = 3.5 respectively. Interestingly, the second of these two 
peaks does not appear at n = 4, which is typically expected 
for a classic de Vaucouleurs profile. The relative strength 
of these two peaks shifts with increasing wavelength, with 
the stronger disk-dominated peak at n = 1 giving way to 
the spheroid-dominated n = 3.5 peak at wavelengths longer 
than the i band. This is believed to be an indicator of the 
shifting in observed stellar population with wavelength, how- 
ever, see Section [5. 3| for further discussion. In addition, the 
centroid of the n = 1 peak at wavelengths longer than the 
i band appears to move towards higher Sersic index values, 
merging into an elongated shoulder of the relatively stable 
n = 3.5 peak. In the K band, the first peak appears to have 
a mean centred on n ~ 1.5. This should not be surprising, 
as optical bands are more likely to probe the young stel- 
lar populations in the disks of galaxies whereas the longer 
wavelengths pick out the older stellar populations within the 
core regions of a spiral galaxy or in elliptical galaxies. Dust 
may also be an issue at shorter wavelengths, blocking light 
from the core regions of galaxies and therefore biasing re- 
covered Sersic indices towards lower values (for example, see 
Pastrav et al.|20T2 l. It is important to note that these data 
are derived from the same r band selected sample of galax- 
ies observed in different wavelengths, and so these relative 
shifts in peak positions represent real variances in observed 
stellar populations, highlighting a wavelength dependence 
on structural measurements. 

The additional peaks at n ~ 0.05 and n = 20 represent 
failed fits. For these galaxies, the fitting procedure drifted 
into an unrealistic parameter space during the downhill min- 
imisation routine employed by GALFIT. Despite attempts 
to force a better fit from the data within GALFITPIPE, the 
fits to the images of these objects remain corrupted, and are 
not appropriate for further analysis. Bad fits occur for many 
reasons. Typical reasons are over-dense regions introducing 
too many free-parameters into the minimisation routine, or 
bad sky subtraction for that region. The corrupt peak values 
of n ~ 0.05 and n = 20 arise due to constraints placed by 
the fitting code GALFIT, and unchanged for the purposes 
of this study. The upper peak is a hard limit, with galaxies 
unable to obtain a Sersic index beyond this value. The lower 
peak is a result of a consistency check within the GALFIT 
code that attempts to force a fit at n > 0.05, hence lead- 
ing to a small distribution around this value. These errors 
that caused these additional peaks are also the cause of those 
found in the ellipticity distribution, discussed further below. 
The density of objects within these failed regions scales with 
wavelength, with the higher-quality bands exhibiting fewer 
cases than poorer quality bands such as the u band. A con- 
servative estimate (n < 0.07 or n > 19) places 1.1% (1,456) 
of r band galaxies within these extremely non-physical re- 
gions, rising to 9.1% (12,630) in the worst affected u band. 



Distributions of recovered effective radii (along the 
semi-major axis), converted to kiloparsecs, are also shown. 
Density profiles at all wavelengths appear relatively smooth, 
approximating a skewed Gaussian distribution. Note that 
these distributions exhibit no additional peaks as observed 
in the Sersic index and ellipticity plots. When regarding the 
median values of these distributions, represented in Figure 
|15| by red dashed lines, we note that the median effective 
radius of a galaxy ranges from 5.5 kpc in the u band to 
3.5 kpc in the K band. This marked decrease in physical 
size with observed wavelength is again as expected if one 
expects the longer wavelengths to probe core stellar popu- 
lations in the bulges of galaxies, whilst shorter wavelengths 
probe recently formed populations in the disks of galaxies 
(i.e., inside-out growth, sec Trujillo & Pohlen (2005)). The 
transition wavelength appears to be the Y band, with a me- 
dian size of ~ 4.5 kpc in the z band, and ~ 3.5 kpc in the 
J band, highlighting once more the wavelength dependence 
on structural measurements and the importance of the right 
choice of wavelength when comparing galaxy samples. In- 
deed, there appears to be little size variation at wavelengths 
longer than the Y band. Clearly, care must be taken when 
comparing the sizes of galaxies observed at different rest 
wavelengths. 

Ellipticity (1 — ^) measurements remain relatively con- 
sistent across all bands, peaking in the range 0.25 < e < 
0.35, and displaying additional peaks at e = and e = 0.95. 
Bands g-K, excepting the r and i bands, show very sim- 
ilar distributions, with a consistent median value of e ~ 0.4 
and a modal value of e ~ 0.35. The higher quality r and i 
bands appear to have, on average, more circular recovered 
ellipticities, with median and modal values of e ~ 0.35 and 
e ~ 0.25 respectively, however the shift is minimal. The 
lower quality u band data leads the ellipticity measurements 
in that band to be biased towards higher values, caused by 
the fitting routine being more susceptible to background 
noise fluctuations and random noise in the frame. Ellipti- 
city measurements presented here are global ellipticities, and 
there will be variability in ellipticity with increasing radius 
from the core on a per-galaxy basis caused by additional 
factors, for example, the effect of seeing or the presence of 
a bar. The additional peaks directly correspond to those 
already discussed previously, and represent failed fits. The 
values of e = and e = 0.95 correspond to internal GALFIT 
boundaries. 



4-3.3 Sersic Magnitudes 

It is not known exactly how the light profile of a galaxy 
behaves at large radii away from the core regions. The ex- 
act nature of any profile will undoubtedly be influenced by 
many factors including, but not limited to, recent merger 
history, star formation rate, gas and dust content and local 
environment. Magnitudes within isophotal radii, not surpris- 
ingly, systematically underestimate the total galaxy light, in 



particular, relative to the Sersic magnitude (e.g., Caon et al. 
1990[|1993|. [Graham fc Driver] (|2005|) show for example that 



Kron magnitudes may underestimate the total galaxy flux 
by as much as ~ 55% dependent upon choice of the mul- 
tiple of Kron radii chosen to integrate out to and the profile 
shape of the galaxy. The comparative value for Petrosian 
magnitudes is considerably worse, underestimating flux by 
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Figure 15. Global results from the SIGMA common coverage sample for all nine bands. Each column represents (from left to right) 
the average model galaxy based on median values for Sersic index, half-light radius and ellipticity and the distributions for all recovered 
Sersic indices, half-light radii (converted into kpc) and ellipticities. The y-axis for each distribution shows the probability density function 
convolved with a rectangular top-hat kernel with standard deviations of 0.05, 0.05 and 0.02 for index, size and ellipticity respectively. 
Median values for each distribution are represented by a red dashed line and are used in creating the average model galaxy in the 
left-hand column. The r band distributions are shown in grey for reference. 
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as much as ~ 95% in the extreme case of a high-Sersic index 
object integrating out to thrice the Petrosian radius. In ad- 
dition to these considerations, below /is = 27 mag arcsec"^ 
environmental effects begin to play a large role in profile 
shape determination. 

In contrast to traditional aperture methods, studies 
have repeatedly shown the strength of Sersic profiling for 
the majority of elliptical galaxies (e.g.,|Caon et al.| (|1993|); 



Graham & Guzman (2003 1; Trujillo et al. (2004 1; Ferrarese 



et al.| ( |2006| ). |Tal fc van Dokkum] ( |2011| support this view- 
point, showing the light profiles of massive ellipticals are 
well described by a single Sersic component out to ~ 8 re, 
with evidence for additional flux beyond these radii possibly 
related to unresolved intra-group light. With regards to disk 
systems, Bland-Hawthorn et al. (,2005 ) use one of the deep- 



est imaging studies of spiral galaxy NGC 300 to show that 
an exponential profile (n = 1) is a good descriptor of its light 
profile out to ~ 17 re. From a sample of 90 face-on late- type 



galaxies, Pohlen & Trujillo ( 2006 1 confirm the accuracy of 
Sersic profiling down to = 27 mag arcsec"^, and suggest 
up to 10% of their sample show evidence for a deviation 
from a standard n = 1 Sersic fit (Type I), instead showing a 
broken exponential profile. These breaks appear in the form 
of either a downbending (Type II; steeper flux drop-off) or 
upbending (Type III; shallower flux drop-off) with increas- 
ing radii. Importantly, this study also suggests this observed 
feature is independent of local environment. 

It is clear that opinion is divided amongst the com- 
munity as to how a galaxy behaves below the typical lim- 
iting isophote, particularly so in the case of a disk galaxy. 
Each of these studies does however suggest a more complex 
structure at large radii than a Sersic profile extrapolated out 
to infinity would imply. In order to account for the lack of 
reliable profile information at large radii, Sersic magnitudes 
require some form of profile truncation so as not to extrapol- 
ate flux into regimes of which we know little. Two schools 
of thought exist in terms of appropriate truncation meth- 
ods, extrapolating flux down to a fixed surface brightness 
limit or integrating under the profile out to a fixed multiple 
of the half-light radii. A constant surface brightness limit is 
more closely related to galaxy gas content, and so has phys- 
ical meaning. However, this method introduces a redshift 
dependence on truncated flux, causing different fractions of 
light to be missed at different redshifts. Truncating at a 
given multiple of the effective radius assumes that the ef- 
fective radius is well understood prior to truncation, which 
owing to the inter-dependency between output Sersic para- 
meters, is not always evident. It does not display any redshift 
dependence however, and is trivial to subsequently recorrect 
if desired. Corrections are typically minor for most galaxies, 
becoming most acute in high-index systems (see Figure |2]). 

A sufficiently large truncation radius must be adopted 
to provide a close estimate of total flux without extrapol- 
ating too deep into the region of uncertainty. SDSS model 
magnitudes employ a smooth truncation at 3re down to zero 
flux at 4re for exponential (n = 1) profiles and 7re down to 
zero flixx at 8re for de Vaucouleurs (n = 4) profiles. We ad- 
opt a sharp truncation radius of lOre for all Sersic indices, 
which corresponds to an isophotal detection limit of jj,,. ~ 30 
mag arcsec"^, the limit to which galaxy profiles have been 
studied. Figure [16] shows the magnitude offset between the 
Sersic profile integrated to infinity and that truncated at 3, 
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Figure 16. Magnitude offset between the Sersic profile integ- 
rated out to infinity and that truncated at a given multiple of 
the effective radius. SDSS model magnitudes (forcing either ex- 
ponential {n = 1) or de Vaucouleurs (n = 4) profile fits, dotted 
vertical lines) employ smooth tapering truncation radii, represen- 
ted by the shaded red and green areas. SIGMA Sersic magnitudes 
within GAMA adopt a sharp lOre truncation radius, blue line. 



7 and 10 re as red, green and blue lines respectively, with 
shaded areas representing the SDSS tapered limits. A 10 
truncation gives a negligible magnitude offset for n = 1, and 
an offset of Am ~ —0.04 for n = 4, with larger corrections 
for higher Sersic indices. Figure [2] middle panel, shows the 
flux contained within 10 re (dashed vertical line) for various 
values of n. Given a 10 re truncation, ~ 100% of the pre- 
truncation flux is retained for n = 1, reducing to 96.1% for 
n = 4. Sersic magnitudes truncated at 10 re for each galaxy 
processed by SIGMA are adopted as the standard Sersic 
magnitude system throughout the remainder of this paper, 
however, full (i.e., integrated out to infinity) Sersic mag- 
nitudes are provided alongside truncated Sersic magnitudes 
in SersicCatvOl for reference. 

Figure [it] shows the offsets between various magnitude 
systems discussed previously as a function of Sersic index. 
When comparing Sersic magnitudes integrated to infinity 
against SDSS Petrosian magnitudes we see the two systems 
are in good agreement until Ur ~ 2, beyond which the mag- 
nitude offset relation begins to turn-off from the Am = 
relation. This trend argues that Sersic magnitudes are re- 
covering an additional ~ 0.4 magnitudes for an n,. = 8 ob- 
ject which would otherwise have been missed by traditional 
photometric methods. However, for the reasons previously 
discussed, this value should be taken as a rough estimate 
of the true amount of flux missed for an object of a given 
Sersic index. Truncating the Sersic index at 10 reduces 
the scale of this turn-off, as expected, keeping the two mag- 
nitude systems in agreement out to Ur ~ 3, however still 
providing some measure of turn-off beyond this point. We 
would expect SDSS Petrosian magnitudes, or indeed any 
aperture-based photometry, to underestimate total flux for 
objects with large wings, and so this result is not surpris- 
ing and a good indication that truncated Sersic magnitudes 
are performing as expected. The final panel in Figure [17] 
compares truncated Sersic magnitudes against SDSS model 
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magnitudes. The SDSS force fit either an exponential or de 
Vaucouleurs profile fit (marked on the figure by dashed grey 
lines) depending upon which profile an individual galaxy 
most approximates. We see clearly the inadequacy of model 
magnitudes when a more comprehensive Sersic magnitude 
is available, with the population of galaxies split into two 
distinct sub-populations based upon their SDSS forced fit. 
For a galaxy at n = 2 for example, the model magnitude 
for the galaxy may be offset from its correct magnitude by 
as much as Am = ±0.3 magnitudes, with larger offsets ob- 
served for high index galaxies. If one constructs a line of best 
fit for each of these two artificial sub-populations the lines 
pass through Am = and n = 1 or 4 as appropriate, con- 
firming that SDSS and SIGMA agree for exponential and 
de Vaucouleurs type galaxies. As highlighted previously, the 
peak Sersic index for the second sub-population does not lie 
at n = 4 but rather at n ~ 3.5. 



5 VARIATIONS IN STRUCTURAL 

PARAMETERS WITH WAVELENGTH 

5.1 Magnitude Comparisons 

The observed nature and form of a galaxy varies de- 
pendent upon the wavelength at which the observation is 
taken. These variations refiect physical mechanisms occur- 
ring within the galaxy, including but not limited to; dust at- 
tenuation and; intrinsic gradients in stellar population, age 
and/or metallicity (e.g., [Block et alT]|1999j |. Below, we show 
the variance in recovered Sersic parameters with wavelength, 
and discuss how this behaviour is characterised. 

Figure [Ts] compares SDSS ugriz Petrosian photometry 
against truncated Sersic magnitudes as recommended in 
Section |4.3.3| as a function of SDSS magnitude. Each row 
represents a different band, with the mode and standard- 
deviation for varying magnitude subsets inset into the left- 
hand column sub-plots. Across each band we see a good 
global agreement between SDSS and recovered Sersic pho- 
tometry at all magnitudes, with the variance between the 
two photometric methods increasing towards fainter mag- 
nitudes as expected. The global total spread is larger in the 
lower quality u band than in the higher quality r band, ran- 
ging from cju =0.72 magnitudes in the former and Or = 0.21 
magnitudes in the latter. This trend should not be surpris- 
ing, as lower quality data presents a unique challenge in re- 
covering 'correct' structural parameters, with larger errors 
expected between different photometric systems for fainter 
galaxies. In all cases, the peak modal values are typically 
less than Am = 0.03 magnitudes, re-enforcing the notion 
of good photometric agreement between these two different 
methods. This is also in good agreement with the offsets pre- 



(u - r)„ 



2.5 



viously laid out in Hill et al. ( 2011 1. The data points are col 



cured according to the recovered Sersic index, and highlights 
the importance of Sersic modelling in recovering accurate 
structural parameters. At all wavelengths, the largest offset 
between SDSS Petrosian and Sersic magnitudes is observed 
in those well-resolved galaxies whose Sersic indices are large 
(n > 4; red data points). These high-index systems are typi- 
fied by being highly centrally concentrated with large exten- 
ded wings, with the fiux in the wings of these galaxies most 
likely to be missed by traditional photometric methods such 
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Hp 

Figure 17. A series of plots displaying offsets between various 
magnitude systems as a function of r-band Sersic index, witli tlie 
data points coloured according to their u-r rest colour, as shown. 
Contours range from the lO"' percentile to the 90"* percentile 
in 10% steps, (top) The Sersic profile integrated out to infinity 
minus the SDSS Petrosian magnitude and; (middle) the Sersic 
profile truncated at lOre minus the SDSS Petrosian magnitude. 
These figures show how Sersic profiling is able to recover flux 
in the wings of galaxies that would otherwise be missed by tradi- 
tional aperture based methods, such as Petrosian apertures, (bot- 
tom) The Sersic profile truncated at lOrg minus the SDSS model 
magnitude. SDSS force fit either an exponential (n = 1) or de 
Vaucouleurs (n = 4) profile fit to attain their model magnitudes. 
This figure shows how model magnitudes provide an inaccurate 
measure of flux for a galaxy whose Sersic index differs signiflc- 
antly from either n = 1 or n = 4. Vertical dashed grey lines at 
exponential (n = 1) and de Vaucouleurs (n = 4) Sersic indices 
are added for reference. 



as Petrosian aperture photometry. By extrapolating the fit- 
ted light profile, Sersic photometry is able to recover this 
missing fiux and provide a more accurate measure of 'total' 
magnitude. Unresolved compact high-index galaxies agree 
well with the Petrosian aperture method. As an example 
of fiux recovery, Sersic photometry for a relatively bright 
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(r = 18) high-indej|^ galaxy observed in the r band may 
recover as much as Ar = 0.5 magnitudes. This missing flux 
recovery may be seen more clearly in the turn-off shown in 
Figure 17 (middle panel), with the scale of the Petrosian 



correction a direct function of Sersic index. In contrast, in- 
termediate and low Sersic index galaxies (n < 4; blue and 
green data points) agree much more closely with SDSS Pet- 
rosian photometry. It is interesting to note that high-index 
galaxies are also those whose Sersic magnitudes have been 
truncated by the largest amount, and so one must take into 
consideration the arguments laid out in Section [4. 3. 3| when 
analysing these systems. 

Figure [19] shows a comparison between GAMA AUTO 
and Sersic magnitudes as a function of GAMA AUTO mag- 
nitude. It should be noted that the GAMA AUTO photo- 
metry presented here is version 2 data and different to the 



version 1 data presented in Hill et al. (2011 1. Briefly, version 



2 photometry employs an updated source detection pipeline 
over a larger area, with a small fraction of version 1 input 
frames discarded due to erroneous data in that region (e.g., 
badly focussed frames). The process used in deriving these 
renewed data is similar in approach to that previously em- 
ployed, a full description of which may be found in Liske et 
al. (2011; in prep.). As in Figure 18 there is good agreement 
between the two photometric systems, with the larger mag- 
nitude offsets observed in the resolved high-index systems. 



5.2 Two Distinct Populations 

Figures [18] and [19] show that Sersic index plays an import- 
ant role when considering magnitudes, with higher index 
galaxies typically recovering more missing flux than their 
lower index counterparts when compared to traditional pho- 
tometric methods. In addition to this, there appears to be a 
wavelength dependence on the flux difference between high 
and low index galaxies, which has implications for variations 
in other structural parameters with wavelength. In order to 
further analyse this wavelength relationship with structural 
parameters, we deflne two galaxy sub-populations based on 
Sersic index and u — r rest frame colour (AUTO aperture 
deflned, as found in the GAMA catalogue StellarMassesv03 
as described in 



Taylor et al. 2011 1. Figure 20 shows this 



relation for the iC'-band Sersic index, with the data points 
coloured according to stellar mass. The bulk of the galaxies 
appear to lie in two distinct populations, the nature of which 



have most recently been explored in Baldry et al. ( 2006 1 ; 



Driver et al. (20061; Allen et al. (20061; Cameron & Driver 



( 2009 1; Cameron et al. ( 2009 1; Mendez et al. ( 2011 1, amongst 



others, and are typically well described by two overlap- 
ping Gaussians. Blue low-index systems correspond to late- 
type disk-dominated galaxies and red high-index systems to 
early-type spheroid-dominated galaxies. This is well suppor- 
ted by galaxy stellar mass, with the least massive galaxies 
appearing disk-dominated, and the most massive appear- 
ing spheroid-dominated, as expected. The faintest types of 



We note here that by high-index we're referring to n ~ 8 
objects, however, |Caon et al. ( 19931 show that for a deep dataset 
and employing good sky-subtraction methods it is possible to find 
galaxies whose central concentrations are of the order n ~ 15. 
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Figure 18. Comparison between Sersic magnitudes truncated at 
lOre and SDSS Petrosian magnitudes for the five SDSS band- 
passes as a function of SDSS Petrosian magnitude, with the data 
points coloured according to their Sersic index in that band (left 
column) . Vertical lines define subsets at magnitudes brighter than 
those values, with corresponding statistics for mode and standard- 
deviation inset into the figure. Density plots (right column) show 
the relative density of objects in Am-space for each of the afore- 
mentioned subsets. 

galaxy, namely dwarf systems (dE, dSO), are not represen- 
ted in our common coverage sample, and so these two peaks 
do not relate to those morphological classes. We used the 
positions of peak object density for each sub-population to 
deflne a dividing line between them, speciflcally, the line 
which lies perpendicular to one connecting the two peaks in 
object density, bisecting it at the point of lowest object dens- 
ity along the connecting line. The equation of the dividing 
line is given by: 



(n- 



-0.59 log riA' + 2.07 



(9) 



In order to avoid any potential misclassiflcations due 
to the effects of dust attenuation, our longest wavelength 
K-h&nA data was chosen as a measure of central concen- 
tration. Sersic indices recovered at shorter wavelengths re- 
turn a steeper dividing line, with the gradient only becoming 
stable at wavelengths longer than the 2 band. This effect is 
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Figure 19. Comparison between Sersic magnitudes truncated at 
lOre and GAMA AUTO magnitudes for all wavelengths as a func- 
tion of GAMA AUTO magnitude, with the data points coloured 
according to their Sersic index in that band (left column). Ver- 
tical lines define subsets at magnitudes brighter than those values, 
with corresponding statistics for mode and standard-deviation in- 
set into the figure. Density plots (right column) show the relative 
density of objects in Am-space for each of the aforementioned 
subsets. 



characteristic of the effects of dust, and shall be explored 
in more depth in Section |5.3| Interestingly, |Mendez et al.| 
( 2011 1 show that the choice of bands used to quantify colour 
is less important, and so a standard u — r colour definition 
is employed in Equation [9] for comparison with much of the 
current literature. Objects bluer than this dividing line will 
be referred to as disk-dominated late- type galaxies (LTGs), 
whereas objects redder than this line will be referred to as 
spheroid-dominated early-type galaxies (ETGs) throughout 
the remainder of this paper. It is well known that two 2D 
Gaussians are able to aptly describe these two populations. 
It follows therefore that a harsh cut of this nature will no 
doubt introduce a small fraction of cross-contaminants for 
galaxies occupying a parameter space in close proximity to 
this dividing boundary, namely, those galaxies that lie in 
the wings of the opposing Gaussian function. The amount 
of contamination will be small however, with the overall 
trends entirely sufficient for analysing global trends within 
each sub-population. Improvements to the nature of auto- 
matic morphological classification based on global structural 
measurements exhibited in this paper will be the focus of fu- 
ture studies presented in Kelvin et al. (2011; in prep.). 

5.3 Variations with Wavelength 

The key galaxy measurements produced by SIGMA, in 
addition to improved object centring accuracy are: posi- 
tion angle; ellipticity; Sersic magnitude; Sersic index; and 
half-light radius. Understanding how each of these para- 
meters varies with wavelength is crucial to remove biases 
when comparing measurements made in different band- 
passes. Wavelength bias may also represent real physical bias 
caused by dust attenuation and stellar population gradients. 

Recovered position angle should show little variance 
with wavelength, instead varying mainly as a function of 
input data quality. In line with the cosmological principle, 
recovered position angle should merely be a random quantity 
assuming no detector bias, although for small area samples it 
may be coupled with filamentary structure. On a per-galaxy 
basis, one might expect minor variations with wavelength to 
occur in the presence of stellar population gradients in tran- 
sient local features such as star-forming regions and bars, 
with different bands being more sensitive to different stel- 
lar populations that trace distinct structural components. 
However, for the SIGMA common sample of 138, 269 galax- 
ies one does indeed find no noticeable trend with wavelength 
for recovered position angle. 

Recovered ellipticity remains relatively stable at all 
wavelengths, instead varying primarily as a function of the 
quality of the input data, as shown in Figure [15] The highest 
quality r and i bands typically return the most circular 
galaxy models, whereas the lowest quality u and z bands re- 
turn more elongated profiles across the same galaxy sample. 
This is as expected as one reduces the signal-to-noise of the 
data, with the fitting routine becoming increasingly sens- 
itive to nearby background noise, however, further studies 
and deeper data are required in order to comment further 
on this effect. 

Finally, recovered Sersic magnitude is expected to vary 
as a function of wavelength as per each galaxy's individual 
SED, the theory of which is well understood and will not 
be discussed further. This leaves the apparent central con- 
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Figure 20. K band Sersic index versus u — r rest frame colour, 
with the data points coloured according to their galaxy stellar 
mass estimates, as shown. Contours range from the 10*'' to the 
90"" percentile in steps of 10%. The two highest peaks in object 
density, corresponding to two distinct galaxy populations (late- 
type and early- type for low-index and high- index respectively), 
are represented by filled black triangles. The diagonal line lies 
perpendicular to the line connecting these two peaks, and bi- 
sects it at the point of lowest object density along the connect- 
ing line, marked on the figure with a plus sign. This dividing 
line defines two sub-samples, which for the most massive galax- 
ies, relate to disk-dominated systems below the line (LTGs) and 
spheroid-dominated systems above (ETGs), the equation of which 
is inset into the Figure. 

centration (Sersic index; n) and size (half-light radius; Ve) 
of each galaxy as a function of wavelength to be discussed. 



5.3.1 Sersic index with wavelength 

Figure [21] shows the recovered Sersic indices for each galaxy 
in the SIGMA matched coverage sample at their rest-frame 
wavelength, coloured according to their population classific- 
ation as described in Section [K2| Considering the population 
definitions are based on K band Sersic indices it is reassuring 
to note that the spheroidal population primarily retain their 
high Sersic index values across all wavelengths, and similarly 
for the disk population. This indicates a significant level of 
consistency in recovered parameters with wavelength, i.e.; a 
galaxy that appears disk-like in the g band is likely to ap- 
pear disk-like in the H band, for example. 3(j clipped mean 
Sersic indices are shown for each population in each band, 
represented by large filled circles coloured as appropriate. 
Polynomial fits to these mean data points, excluding u band 
values due to their lower quality imaging data, reveal general 



trends in Sersic index with wavelength. The best fit Sersic 
index for the disk-dominated population is given by: 

log rzdisk = -0.715 log' A.est + 4.462 log A.e^t - 6.801 (10) 

and similarly for the spheroid-dominated population: 

logn.ph = -0.2101og' A,est + 1.3941og A,,,t - 1.753 (11) 

where Arcst is the rest-frame wavelength of the observation 
of the galaxy. It is important to remind the reader to be 
mindful of our sample selection when considering these re- 
lations. Note that we have adopted log-quadratic relations 
for Equations [To] and |11| Whilst the spheroid-dominated 
population may not appear to require a quadratic fit, the 
disk-dominated population most-likely does. For this reason, 
the functional form of both equations has been kept the 
same. The linear relation for the disk-dominated population 
is given by: 

log Wdisk = 0.267 log Arcst - 0.676 (12) 

and the spheroid-dominated population is given by: 

lognsph = 0.170 log Arest +0.024 (13) 

These linear relations are provided for reference only and 
are not used in any subsequent calculations, with the log- 
quadratic forms instead being the preferred descriptors of 
the two populations. 

We find that the spheroid population Sersic indices re- 
main relatively stable at all wavelengths, exhibiting slightly 
lower Sersic indices at shorter wavelengths and becoming 
essentially stable at wavelengths longer than the z/Y in- 
terface. Mean ETG Sersic indices range from rig = 2.79 to 
riK = 3.63 from g through to K, an increase of 0.11 dex, 
equivalent to 30%. This increase is consistent with the 23% 
increase reported in |La Barbera et ah] ( |2010| over a similar 
wavelength range. However, whilst the fractional increase is 
comparable, the absolute values are not; |La Barbera et aT] 
find on average Sersic indices n ~ 2 — 3 larger than those re- 
ported here. Whilst it is unclear what causes this difference, 
a potential difference in sample definitions may be import- 
ant. That study defines ETGs based on a number of SDSS 
parameters including fracDeVr; a parameter that describes 
how well the global light profile of the galaxy is fit by a de 
Vaucouleurs profile. A cut of this nature is somewhat ana- 
logous to making a Sersic index cut alone which, as can be 
seen in Figure p^anda gain in Figure [20l would introduce an 
element of contamination from the LTG population. If a rel- 
atively large number of the ETG sample in |La Barbera et al.| 
are in fact bulge-dominated systems with a weak underlying 
disk then one might expect the Sersic indices of their bulges 
to differ somewhat from a traditional de Vaucouleurs profile. 
A Sersic index of n ~ 6 is the value found in Simar ~et al.l 
( |2011[ ) for bulge-|-disk systems with a well-defined bulge, in 
good agreement with the offset found here. 

The apparent stability found in Sersic index with 
wavelength is as expected for relatively dust-free single- 
component early-type systems, and is interesting to re- 
confirm empirically. Since the spheroid-dominated popu- 
lation is likely to include a small fraction of misclassi- 
fied galaxies, as previously discussed, due to the nature 
of the harsh cut presented in Section |5.2[ a small gradi- 
ent with wavelength should not be surprising. Recent work 
by Rowlands et al. (2011; submitted) suggest that as many 
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as 5.5% of early- type galaxies contain significant fractions 
of previously unaccounted-for dust, introducing an addi- 
tional secondary deviation in recovered Sersic indices with 
wavelength. Dust in a galaxy is typically centrally concen- 
trated, and so has the effect of masking stellar light eman- 
ating from the core regions of a galaxy. Since galaxy fitting 
algorithms such as GALFIT apply larger weighting to higher 
signal-to-noise regions, minor deviations at small radii have 
the potential to drastically affect the recovered structural 
parameters, including the Sersic index. Therefore, the addi- 
tion of dust to the core region of a galaxy would subdue the 
cuspiness of the galaxy and bias the model towards a lower 
Sersic index. 

The disk population exhibits a larger change in Ser- 
sic index variation with wavelength than that observed for 
the spheroid population. The recovered mean disk Sersic in- 
dices range from Ug = 0.92 to tik = 1.40 from g through 
to K, an increase of 0.18 dex, equivalent to 52%. As with 
the spheroid population, disk Sersic indices also become in- 
creasingly stable at wavelengths longer than the z/Y inter- 
face. Since we typically expect disk-dominated systems to be 
dustier than their early- type counterparts, owing to the pre- 
valence of ongoing star- formation in many of these galaxies, 
then a significant variation in Sersic index with wavelength 
should be expected as a consequence of the arguments pre- 
viously laid out. Since the disk Sersic index appears stable 
beyond the z/Y interface, we can conclude that the effect 
of dust in these regimes is minimal, and therefore if 'in- 
trinsic' disk Sersic indices are required, one should look to 
the longest wavelength data available, typically longwards 
of the z band. In addition to the effects of dust attenuation, 
we may also consider stellar population gradients. Since this 
sample is not a pure-disk sample, and instead contains a 
host of disk-dominated systems, a fraction of galaxies in 
the disk-dominated population will no doubt contain ad- 
ditional structures such as a bulge and/or a bar. Bulges 
tend to contain older, redder stars of a higher metallicity 
than the younger, bluer stars found in the disks of galaxies. 
Shorter wavelengths are more sensitive to the blue popula- 
tion found in the disk whereas longer wavelengths become 
increasingly sensitive to the red population. Therefore, any 
real colour gradients that exist in a galaxy, which are indic- 
ative of metallicity and age gradients in the underlying stel- 
lar population distribution, would also lead to a change in 
the measured Sersic index, dependent upon the wavelength 
at which that galaxy was observed. A short wavelength is 
therefore more likely to probe the disk stellar population 
than a longer wavelength. It is unclear whether the effect 
of dust attenuation or stellar population gradients are the 
dominant factor in determining the variation in Sersic in- 
dex with wavelength, with a combination of both likely to 
contribute globally. 

We note that the disk-dominated and spheroid- 
dominated populations, once stabilised, tend towards 
n-disk — >■ 1.4 and risph — >■ 3.6 respectively. These values dif- 
fer from the Sersic indices typically used to describe late 
and early-type systems (excluding dwarf galaxies, for which 
there is a Magnitude-Sersic index relation), namely niato = 1 
and nearly = 4 respectively (represented in Figure 21 by hori- 



zontal grey lines) . This may indicate morphological contam- 
ination between populations as previously discussed, with 



some galaxies exhibiting bulge-to-disk ratios away from val- 
ues of either zero or unity. 



5.3.2 Half-light radius with wavelength 

Figure[22]displays the recovered half-light radii as a function 
of their rest-frame wavelength. 3a clipped means are repres- 
ented by solid red and blue circles for the disk-dominated 
and spheroid-dominated systems respectively, with linear 
fits to these data shown. The best fit linear relation de- 
scribing the half-light radii in physical units (kpc) for the 
disk-dominated population is given by: 



logrc.disk = — 0.1891og Arcst + 1.176 
and for spheroid-dominated systems: 
logre,sph = -0.304 log Arcst + 1.506 



(14) 



(15) 



where Arcst is the observed rest-frame wavelength for the 
galaxy. Again, it is important to remind the reader to be 
mindful of our sample selection when considering these re- 
lations. 

Using these relations, we observe significant variation in 
the recovered sizes of galaxies as a function of wavelength. 
The disk population mean half-light radii range from r^.g = 
4.84 kpc to re,K = 3.62 kpc from g through to K, a decrease 
in size of 0.13 dex, equivalent to a drop of 25%. The spher- 
oid population exhibits a larger spread from r^.g = 5.27 kpc 
to re,K = 3.29 kpc over the same wavelength range, a de- 
crease in size of 0.20 dex, equivalent to a drop of 38%. This 
variation in the spheroid population size is in good agree- 



ment with studies by La Barbera et al. (2010 1 ; Ko & Im 



(20051, reporting decreases of 29% and 399(|^°| respectively 



over a similar wavelength range. Explanations for the vari- 
ation in recovered size with wavelength include dust attenu- 
ation at shorter wavelengths or metallicity gradients within 
the galaxy. The effects of dust on a galaxy light profile have 
previously been discussed. Obscuring the central region of a 
galaxy would shift the balance of total flux towards larger 
radii, artificially increasing the half-light radius. This effect 
is well understood for late- type systems, and indeed has been 



predicted in the literature notably by Evans ( 1994 1; Cunow 



(2001 1, and more recently by MoUenhoff et al. (2006 1; Gra- 



ham &: Worley| (20081. Data from these studies are added 
into Figure |22| for reference (as indicated), normalised to 
the disk half-light radius predicted at Aobs = 900 nm. Where 
possible, we adopt or infer a dust face-on optical depth of 
3.8 (central face-on B-band opacity) as recommended 
Driver et al. ( 2007 1 . This value is close to that found in de- 



tailed modelling of nearby galaxies by Popescu et al. 
(|200T| 



Misiriotis et al. 



Popescu et al. (2011 1; Hermelo et al 



(20001 



(in prep.). We also adopt an average inclination galaxy of 
cos^ = 0.5. |Evans| analyses face-on galaxies alone, whereas 



Cunow 



only produces detailed models for = 3.0 systems, 
and so care should be taken when comparing these data. 
In all cases we see an excellent agreement between the ob- 
served size- wavelength variation and that predicted by these 
dust model simulations, notably so when compared with the 



A linear extrapolation of the trends reported in |Ko fc fm| 
l |2005l l were applied in order to convert V-K offsets into the g- 
K wavelength range used here. 
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Figure 21. Recovered Sersic index shown as a function of rest-frame wavelength in log-log space, coloured according to the population 
definitions described in Section |5.2| Blue data-points correspond to disk-dominated galaxies whereas red-data points correspond to 
spheroid-dominated galaxies. Large red and blue circles show the 3(j-clipped mean Sersic indices for each respective population in each 
band, positioned at the median-redshift rest-frame wavelength for that population. Polynomial fits to these mean Sersic indices are shown 
for both populations, the equations of which are inset into the figure. Owing to its lower quality imaging data, we exclude the u-band 
data in the calculation of these lines. Vertical lines show the Ic spread in the data. For reference, grey horizontal lines at n = 1 and 
n = 4, equating to exponential and de Vaucouleurs profiles respectively, are added. 



work of |Mollenhoff et al.| as they employ dust models that 
can account for both dust attenuation and emission. As a 
caveat, we note that the population shown in Figures |21| 
and [22] does not constitute a volume-limited sample, and so 
conclusions must remain tentative until further studies can 
confirm this relation. However, our initial interpretation is 
that dust models more than adequately account for the ap- 
parent size-wavelength relation in late-type disk-dominated 
galaxies. 

It is interesting to note that whilst the spheroid pop- 
ulation shows less variation in central concentration, i.e., 
Sersic index, than the disk population, it exhibits a larger 
size variation with wavelength. Dust is not expected to be 
a dominant factor in the attenuation of hght within these 
systems (although see the earlier discussion regards early- 
type dust fractions). However, higher optical depth values 
than that recommended in Driver et al. (2007l would have 



the effect of skewing the gradient of the dust attenuated 
size-wavelength relation to match the observed distribution. 
In addition to the possibility of age/metallicity gradients 
within spheroids, an alternative explanation for the apparent 
size variation with wavelength in early-type systems relies 
upon the interdependency between recovered Sersic index 
and half-light radius. A change in the Sersic index arising 



due to, e.g. small core dust components, additional unre- 
solved or disturbed structure in the core arising due to re- 
cent environmental interactions, the influence of an AGN 
or uncertainty in the PSF may lead to an equivalent cor- 



rective change in half-light radius. MacArthur et al. ( 2003 1 



show that uncertainty on the PSF fuU-width-half-maximim 
by as much as F = 1.5" ± 0.5, a range encompassing most 
of the optical data as shown in Figure |10| would yield an 
equivalent measured size variation of 25% for the worst af- 
fected compact systems. One concern might be that small 
numerical uncertainty in Sersic index would yield artificial 
changes in recovered size. Using the Sersic index ranges for 
both disk and spheroid populations shown in Figure [21} and 
fixing the effective surface brightness and total magnitude 
in the total luminosity variant of Equation 2 in Gr aham fc| 
[Driver (20071, we would expect to see equivalent changes in 
Sersic half-light radii of Arg^disk = 17% and Are_sph = 8% 
for the disk and spheroid populations respectively, i.e., far 
less than that seen here. 



5.3.3 Co-variation of Sersic index and half-light radius 

We have shown how the Sersic index and half-light radius 
for the spheroid and disk populations vary with wavelength. 
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Figure 22. Recovered half-light radii in kpc as a function of rest-frame wavelength. Red and blue circles show the Scr-clipped mean half- 
light radii for spheroid-dominated and disk-dominated galaxies respectively in each band, positioned at the median-redshift rest-frame 
wavelength for that population. Linear fits to these mean half-light radii are shown for both populations, the equations of which are 
inset into the figure. Owing to its lower quality imaging data, we exclude the «-band data in the calculation of these lines. Overlaid are 
data from several authors who predict an increase in the measured half-light radii in late-type systems due to the effects of dust. Further 
details are available in the text. 



however, one must not consider these variations in isolation. 
All of the output model parameters have a combined ef- 
fect on the final light profile of a galaxy. Several of these 
parameters including: Sersic index; half-light radius; total 
magnitude; and the background sky display certain levels 



of inter-dependence (e.g., Caon et al. 1993 Graham et al 



1996 1. For example, an over-estimation of the sky level would 
lead to an under-estimation in the total magnitude of that 
galaxy, and consequently Sersic index and half-light radius 
also. In the case of sky however, the signal-to- noise weighting 
employed by GALFIT somewhat ensures against sky offsets 
of a few counts from the true value noticeably adversely 
affecting the fit. Since we would not expect the system- 
atic error in the background sky to show significant trends 
with wavelength, and the variation in total magnitude with 
wavelength is well understood, we exclude them from our 
investigation into the co-variation of structural parameters 
with wavelength. We now consider the combined effect of 
varying the Sersic index and half-light radius in unison, and 
how this impacts on the overall hght profile of a test galaxy 
from u K. 

Using Equations [Tol [TT| [U] and [Ts] we generate estim- 
ates of Sersic indices and half-light radii at equal steps in 
log-wavelength space for both the spheroid and disk popu- 
lations. Using the Sersic relation, and assuming a constant 



total magnitude for both spheroids and disks of mtot = 15, 
we create a series of surface brightness hght profiles from 
u through to K. Figure [23] shows the change in the re- 
covered surface brightness light profiles over the u ^ K 
wavelength range, with the shaded areas representing the 
maximal area swept out by these hght profiles as they vary 
in wavelength. This gives us an indication of how changes in 
recovered structural parameters affect the underlying sur- 
face brightness profile. The hatched region represents the 
worst-case limit at which these hght profiles may be accep- 
ted as containing significant signal above the background 
sky level, as given in Section [4. 2. 4| The vertical dashed line 
represents a 1 pixel distance from the centre. 

Despite the relatively large size variation observed in 
the spheroid population (a decrease of 38% in g — > K) , when 
considered in conjunction with the Sersic index variation (an 
increase of 30% in g — > K) the combined effect amounts to 
a relatively modest impact on the majority of the recovered 
light profile. It appears that as the spheroidal size decreases 
the Sersic index increases at a comparative rate. The most 
noticeable surface brightness variation is found in the cent- 
ral core region, fluctuating by 0.49 magnitudes at the 1 pixel 
boundary. Since a significant fraction of total fiux lies in the 
core regions of high-index systems, it should not be surpris- 
ing that a small variation in Sersic index would produce a 
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relatively large variation in half-light radius. Despite this ef- 
fect, the majority of surface brightness profile out to large 
radii remains relatively stable with wavelength, vastly redu- 
cing the need for more complex mechanisms as previously 
discussed. 

The variation in size for the disk-dominated population 
(a decrease of 25% in (? — >■ K) coupled with a relatively 
large increase in Sersic index (an increase of 52% in g — > K) 
yields a similar effect on the surface brightness profile vari- 
ation as previously described for the spheroid population. 
Surface brightness fluctuates by ~ 0.86 magnitudes at the 1 
pixel boundary, an increase of 75% on the variation in the 
spheroid population. Here it appears that the impact of dust 
attenuation has a particularly distinct effect on the light pro- 
file in disk-dominated galaxies, agreeing well with the the- 
oretical predictions for size variation with dust presented in 
Section I5X2] 

Whilst no single mechanism can be shown to be entirely 
responsible for the relations between Sersic index, half-light 
radius and wavelength observed across the two populations, 
it is clear that the large apparent size fluctuations in the 
spheroid population appear to be initially misleading. Only 
when considering Sersic index in conjunction with half-light 
radius does the true nature of these effects come to the fore. 
The spheroid population, despite exhibiting large changes 
in half-light radius with wavelength, maintains a relatively 
stable surface-brightness profile from u through to K. The 
variation in the disk population with wavelength appears 
well described by current dust models, however, it is most 
likely a combination of dust attenuation, stellar population/- 
metallicity gradients, unresolved secondary features in the 
core region affecting profile fits, and uncertainty on addi- 
tional parameters such as the PSF that affect the underlying 
physics in these systems. Future studies aim to further in- 
form this discussion for a limited sub-sample to be presented 
in Kelvin et al. (2011; in prep.). 
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Figure 23. Surface brigiitness variation from u througii to K for 
tlie early-type splieroid-dominated and late-type disk-dominated 
populations. We generate Sersic indices and half-light radii in 
wavelength steps from u K for each population using the 
trends as described in equations |10[[TT]|14| and |15[ Using these val- 
ues, surface brightness profiles (without PSF convolution) may be 
constructed for each wavelength bin. The shaded regions shown 
here represent the maximal area swept out by these light pro- 
files along the transition from u through to K, and represent how 
much of an effect the reported changes in Sersic index and half- 
light radius have on the overall light profiles. The hatched region 
indicates the brightest limit at which light profiles may be trusted 
(Mlim.if = 22.07 mag/arcsec~^), and the vertical dashed line rep- 
resents 1 pixel in distance from the centre. Profiles are produced 
assuming a constant total magnitude of mtot = 15 for both the 
spheroid and disk populations. 



6 CONCLUSION 

We have produced high-fidelity automated two-dimensional 
single-Sersic model fits to 167, 600 galaxies selected from the 
GAMA input catalogue. These have been modelled inde- 
pendently across ugrizY J H K using reprocessed SDSS and 
UKIDSS-LAS imaging data. These data have subsequently 
been delivered to the GAMA database in the form of the 
catalogue SersicCatvOl. In order to facilitate the construc- 
tion of this dataset, SIGMA, an extensive multi-processor 
enabled galaxy modelling pipeline, was developed. SIGMA 
is a wrapper and handler of several contemporary astronomy 
software packages, employing adaptive background subtrac- 
tion routines and empirical PSF generation on a per-galaxy 
per-band basis to tailor input data into the galaxy model- 
ling software GALFIT 3. Output results from GALFIT are 
analysed for pre-determined modelling errors such as pos- 
itional migration, extreme model shape and/or size para- 
meters and adverse nearby neighbour flux. Nearby object 
masking is employed as a last resort, with secondary neigh- 
bours being preferentially modelled simultaneously with the 
primary galaxy in the first instance. 

Using this dataset, we have defined a common cover- 
age area across the three GAMA regions that encompasses 



138, 269 galaxies, 82.5% of the full sample. This common 
area contains only those galaxies which have been observed 
in all nine bands, providing a useful basis upon which to 
further explore wavelength trends. We define a Sersic mag- 
nitude system that truncates Sersic magnitudes at 10 r^- 
This ensures that flux is not extrapolated below the typical 
limiting isophote into regions where data quality and quant- 
ity is not sufficient to constrain the form of the galaxy light 
profile. Truncated Sersic magnitudes appear to be a good 
descriptor of global galaxy colours and total galaxy fiux. For 
well-resolved disk-like galaxies (n < 2), traditional aperture- 
based methods are in good agreement with truncated Sersic 
magnitudes. For high centrally-concentrated systems how- 
ever (n > 4), it appears that traditional aperture-based, 
such as Petrosian magnitudes, may miss as much as Am^ = 
0.5 magnitudes from the total fiux budget which is only re- 
covered through Sersic modelling. 

When considering the dataset in n-colour space we find 
galaxies appear to exist in two distinct groups. For the 
most massive systems, we associate these two groups with 
the spheroid-dominated early-type galaxy (ETG) and disk- 
dominated late-type galaxy (LTG) populations. Owing to 
the nature of our input sample selection, these definitions 
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do not extend down to the fainter dwarf population, and so 
subsequent trends will not represent those systems. We use 
the longest wavelength K band Sersic index measurements 
in conjunction with rest-frame u — r colour to define these 
two populations. Using these definitions, we are able to fur- 
ther probe the variations in recovered structural parameters 
with wavelength for each population. 

We find that the Sersic indices of ETGs remain reas- 
onably stable at all wavelengths, increasing by 0.11 dex 
(-1-30%) from g to K and becoming very stable beyond the 
z/Y interface. In contrast to this, we find that LTGs exhibit 
larger variations in Sersic index with wavelength, increasing 
by 0.18 dex (-1-52%) across the same wavelength range. Re- 
covered sizes for both the spheroid and disk systems show 
a significant variation with wavelength, showing a reduction 
in half-light radii of 0.20 dex (-38%) in ETGs and 0.13 dex 
(—25%) in LTGs from g to K. Size variation of this scale 
for disk systems has been well predicted by dust models, 
highlighting the important role dust attenuation plays when 
considering structural variations across a broad wavelength 
range. 

We note that spheroidal systems exhibit a larger size 
variation with wavelength than that found in disk systems. 
Possible physical explanations for this behaviour include low 
levels of unresolved dust or the effects of AGN feedback in 
the core of the galaxy, both of which would affect Sersic 
profiling. Significant amounts of dust, such as an increased 
dust attenuation optical depth parameter r^, may allow 
current dust models to accurately describe the variation in 
half-light radii we find. It is unlikely however that a signi- 
ficant fraction of our spheroid-dominated population con- 
tain sufficient amounts of dust for this to be the case. Large 
stellar population/metallicity gradients present within indi- 
vidual structures of the galaxy would cause galaxies to look 
markedly different in different wavelengths, contributing to 
any concentration- wavelength/size-wavelength variation. In 
addition to these factors, uncertainties on the measured PSF 
and background sky must be considered. 

However, when considering variations in half-light ra- 
dius and Sersic index together with wavelength we find that 
the large fluctuations in spheroidal parameters amount to 
a relatively modest impact on the recovered light profile. 
A comparatively larger effect is noted for the disk systems, 
particularly in the core region, supporting the presence and 
effect of dust attenuation in addition to stellar population/- 
metallicity gradients. At a distance of 1 pixel from the cent- 
ral region, spheroid systems display a variation in surface 
brightness of 0.49 magnitudes from u through to K. In 
disk systems, the comparative figure is 0.86 magnitudes, 
an increase of 75%. This highlights the importance of not 
considering recovered parameters in isolation, as the inter- 
play between them has the possibility of masking underlying 
trends. 

The effects of dust attenuation appear to be the dom- 
inant factor constraining the variations in structural para- 
meters with wavelength, notably so for the disk-dominated 
population. In contrast with this, apparent large structural 
variations in the spheroid-dominated population appear to 
have a relatively minor effect on the underlying surface- 
brightness profile than might have been expected. Future 
studies in Kelvin et al. (2011; in prep.), focussing on a lim- 
ited sub-sample of this dataset, will provide a deeper un- 



derstanding of these structural variations with wavelength, 
enabling us to comment further on the key mechanisms in- 
volved in varying structural parameters with wavelength for 
a host of different morphologies. 
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APPENDIX A: SIGMA INPUT OPTIONS 

On starting SIGMA, a number of input options can be spe- 
cified. Some of these are essential in its use, whereas others 
are designed for testing purposes only. The available input 
options can be found in the help document, reproduced be- 
low. 

Listing 1: SIGMA help lists the available input options that 
may be specified when starting SIGMA. 



SIGMA Ve 



0.9-0 — 23 Jul 2010 



DESCRIPTION 

SIGMA (Structural Investigation of Galaxies 
dimensional fitting code taking inputs from 
producing models using the GALFIT software. 



via Model Analysi 
the GAMA SWarped 



OPTIONS 



make a plot o 
output catalo 



# # 

# ,# 



r of bootstrap 
mple , from lowe 
iple allowed mu 



thod (1=GALFIT, 2=B00TSTRAP) 

ut .fits files (png format) 
tput catalogue , only data 



in GALFIT 
1,2,5,10) 



GAMA ID 
SIGMA ID 
SDSS OBJID 



In this study, we initialised SIGMA using the following 
command: 

sigma — b x — n sigmacat x . csv — p 16 — t 1 

where x represents the band for modelling (ugrizY J H K). 
This command initialises SIGMA on 16 processors, restrict- 
ing the fits to single-component (single-Sersic) only. Running 
SIGMA across each band individually, we produced 9 indi- 



append A to output log files 

A-band (default: r) 

input catalogue [img/csv] (ne 



RA & DEC) 



vidual catalogues for later matching in TOPCAT (Taylor 
l2005|. 



APPENDIX B: INITIAL CONDITIONS 

The galaxy modelling phase has been discussed extensively 
in Section|3] Most of the input parameters fed into GALFIT 
come directly or trivially from the pipeline OBJECTPIPE (Sec- 
tion [3!5]), a module wrapper around Source Extractor. The 
two main exceptions to this are half-light radius and Sersic 
index. 

Half-light radii from Source Extractor are modified be- 
fore being fed into GALFIT. This is to account for the differ- 
ence in radii definitions between the two programs. Source 
Extractor's FLUX RADIUS parameter outputs a circular- 
ised radius which is based on PSF convolved imaging data. 
The format of GALFIT's initial estimate of the half-light 
radius is that along the semi-major axis which is intrinsic 
to the object (i.e. - deconvolved from the PSF). Equation [5] 
converts Source Extractor circularised radii into semi-major 
intrinsic radii appropriate for GALFIT. Figure [Bl] displays 
the before (uncorrected) and after (corrected) Source Ex- 
tractor half-light radii against their output modelled half- 
light radii from the SIGMA full sample, coloured according 
to their predicted morphological type as detailed in Section 
|5.2| Unmodified Source Extractor half-light radii are a poor 
predictor of final modelled GALFIT half-light radii, as ex- 
pected. Large galaxies (re > 4 pixels) tend to have their sizes 
underestimated by Source Extractor by as much as 50%. Fol- 
lowing a turn off at ~ 4 pixels, small galaxies tend to have 
their sizes overestimated by Source Extractor. Once these 
data have been corrected, we find a marked increase in the 
agreement between the two measures, notably so for late- 
type galaxies (blue data points). Data above re > 4 pixels 
has been reduced to minimal scatter about a 1:1 correlation. 
The effect of the turn-off has been significantly mitigated, 
yet not entirely diminished. This indicates the difficulty in 
accurate size estimation of galaxies that only occupy a mat- 
ter of a few pixels. Correcting radii in this manner signific- 
antly reduces the chance of GALFIT finding a local-minima 
during the minimisation phase, and consequently reduces 
the risk of convergence on a non-physical solution. 
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FLUX RADIUS / pixels FLUX RADIUS (modified) / pixels 

Figure Bl. A comparison between Source Extractor half-light radii (FLUX RADIUS) and modelled GALFIT half-light radii (re) in 

the r-band, with data points coloured according to their predicted morphological type as described in Section |5.2| (Left) Uncorrected 
half-light radii from Source Extractor are a poor initial condition for modelling an object in GALFIT, as Source Extractor radii are 
circularised and make no attempt to correct for the effect of seeing. (Right) Using Equation [s] to correct for ellipticity and the PSF 
produces a much better starting value for GALFIT, hence reducing the chance of finding local minima during the minimisation phase. 
The green line represents a 1:1 ratio, for reference. 



There is no obvious proxy for Sersic index in the de- 
fault Source Extractor parameters file. An approximation 
was created based on the trend between the output Sersic 
index and the ratio between the corrected half-light radius to 
the Kron radius for a small test sample of trusted galaxies. 
From this we were able to derive a relation for a predicted 
variable Sersic index: 



where re is the corrected Source Extractor half-light radius 
(Equation^ and rKron is the Source Extractor Kron radius. 
Figure [B2| sliows the density distributions between the vari- 
able Sersic index, rivar, and several other initial conditions 
for a sample of 49, 395 galaxies in the r band. The other in- 
put parameters (size, position angle, ellipticity, magnitude, 
position) are not modified. These results show that the final 
recovered Sersic index is largely independent of its initial 
condition, with the notable exception of a bump in the dis- 
tribution at n = 0.1 for riinitiai = 0.1 and a variable height 
spike of failed objects at n ~ 20. The n = 0.1 bump repres- 
ents galaxies whose initial Sersic index guess is placed too far 
away from its true value, and so fails to successfully migrate 
away from the initial parameter space using the Levenberg- 
Marquart method employed by GALFIT. It appears that 
the minor fluctuations found in the main body of the dis- 



tributions directly correspond to the varying height of the 
n ~ 20 spike. Despite these features, it is clear that the 
initial Sersic index is afforded a great deal of variability in 
order to achieve a successful and consistent fit. The majority 
of distributions presented in Figure [B2] show little variation, 
with similar levels of success and failure. It was therefore felt 
that a simple ninitiai = 2.5 would be an appropriate initial 
condition as it lies in the middle of the expected probability 
space, yet not at the value of either of the bimodal peaks. 
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Figure B2. A plot comparing final modelled Sersic indices for a sample of 49, 395 galaxies in the r band given different initial Sersic 
indices, as shown. From top to bottom, the initial Sersic indices fed into GALFIT are; variable (see Equation |BiJ|; 0.1; 0.5; 1.0; 2.5; 4.0; 
10.0. Underlying the fixed initial Sersic index distributions is the distribution for the variable Sersic index coloured in grey, for reference. 



© 2011 RAS, MNRAS 000, [TpS] 



